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liTRODyCTlcati 

A whicl® RWrtfing i» cootiHet wim tn® growawl follow® 
tho ground profit®* grouM iaii@%®fiii#s® tliorofor® ®et« o® ® 
hmm «tcit*tli»s for tlwt voMcl® througli tha point® of cwmtoot* 
Thm natiir® of tho ground initocod exisltttion i® a®toaiid.i»d toy 
tho ctiaroctoristie® of tbe grouiSKl onw^iiiiis® ^ typo of €ont»«t 
and wshiolo «p®«d* Mmst wtiicl® €q?®rat«« on prepared traolsi 
«acii a® road®* rail®# rtniway® ®to* Thl# ttie®!® d««ls with 
respons®^ of ntilway wagon®# oasoli®® oto*# to the wertie®! 
onwoiMiiMi® of the rail tracii* 

^ rail txmtM oonsist® of fiohplated rail Joint® ma& 
i® siai^rtad on olaaper® oad ballast*/ 2tia slaapars bold tba 
rail in oorraot aligMtnt and treo^iaiid.t*' tba traidiE looS to 
tba ballast ower a wif® area* A rail Joint is tba wa«ic«8t 
part of tsb® tr^ and a atody of tba affaot of tbasa Joliit® ' 
on aartioal responMi and tradte foroos 1® wary ia^rtant 'frofti 
tba point of wlaif of daoign' of aabiola imd tracic* A wary 
oQMooii faatnra of tba J^nt® i® tba bogging of tba rail® ' 

'the 

witb^ fossaga of wOhiola# i#ian bogging anoaad® a oartsin 
Unit# tba Jointa mm oallad low Jointa#: Ao.'baawy nilntanonoa 
I® ragairad for tbana J^nts# waldad Jointa ora noW' bo.ing 
peowidad owor soaaral wald langtlw on bi#i spaod traac# : 
Howwaar# iii^platad Joint® still, oselst'^' to alloNf tbasnal ^ ' 
oiqponiion of tba rail®.# 
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The tmewnness of the trade can he considered to be 
of two types* 

1. Continuous unevenness due to irregularities of the 
rail profile, and 

2m Discrete unevenness due to hogging of the rails at 
the Joints, 

The study of vibration effects on the vehicle due to 
above two types of unevenness are significant in respect of 
i) design of vehicle structure, il) passenger txanfort, 
ill) controllability of the vehicle, and iv) design of the 
tradk. The vibration of the vehicle Induces dynaiaic loads 
on the vehicle BBcid tradk structure* The ground induced loads 
may cause fracture or fatigue failure of structural components. 

Literature Survey 

The response of vehicles to cemtinuous trade unevenness 
has been studied by several investigators^ bnly a few 
iz:^stigatlons are reported in open literature for discrete 
low Joints, Radford (4) has investigated the vertical fore® 
between rail idieel and rail at dipped rail Joint, % attempting 
a deterministic solution, he has shown that the iiapulsive 
force produced' due to passage- of heavi'ly loaded vehicles on 
sudi irregularities increased linearly with the speed. By 
using syMetrical rail Joints, he has also shosn that the 
re<hicticsi in unsprung mass cause® the reduction in contact 
force idii€h is Boost isoportant for the deterioration of the 
Joint, 
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In reference (14)# Fryba has given various methods for 
calcjulating the vibration effects of base when a moving load 
passes over It* 

In reference ( 9 ), Yadav has extensively caltjulated the 
response of different vehicle models over randcro tracSc 
irregularities by taking unevenness as a homogencxis process* 

Analytical Inyestioationa 

The analytical Investigation reported in thesis covers 
three aspects# 1) mathematical model of the vehicle# ii) model 
of track unevenness# iii) solution of governing equations to 
obtain response statistics* 

Single degree of freedom model with linearised spring 
and daiBpers is investigated. It is assumed that vehicles 
moves with constant velocity. Model of track unevenness 
consists of continuous irregularities of the vertical rail 
profile plus unevenness of the Joints. The governing equations 
are solved by using convolution integral approach* 

The parameters investigated in the thesis are# 
i) maximaiB value of vertical displacement# ii) maximum value 
of vexHtical acceleration* and iii) maximum Mid minimum value of 
the contact force. The upper limit to the response as well 
as rejipoiise acceleration is Important for passenger oaefort 
and drivers ability to ccxitrol the vehicle, Ihis net^ a limit 

to tee vehicle velocity and msilies tee dpiMedc response a 

. ' ' ‘ s 

constraint in vehicle under-carriage design. The upper limit to 
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tihe contact force laay be iaportant to the vdbicle control 
system as well as for stress calculations in structural 
ccm^onents of the trac3c and vehicle. 2!he lower value of contact 
force is important for derailment studies. 

A chapterwise summary of the contents of the thesis 
is given below* 

Chapter 2 gives a description of mathematical imadels 
of the vehicle and track unevenness. derivation of the 

equation of motion and their solution by convolution integral 
have been discussed. 

In Chapter 3* response parameters for harmonic 
excitation are determined. The analysis of such a theoretical 
excitation to the v^icle is useful as it will form a basis 
for Fourier analysis considered later in Chapter 4, Since 
response to harmonic excitation is easy to calailate, Fourier 
analysis Is an isportant tool for calculating responses of 
different type of excitations. 

Chapter 4 describes the response of vehicle to the 
periodically occurring low joints. Response parameters 
are determined by using two methods. <a> single pulse 
excitationi in this analysis response of the vehicle for sin^e 
joint excitation is calculated. Some suitsfele shape of the 
pulse Ijoint) is assumed* Then response for series of pulses 
are obtained by superimposing response of eaCh pulse, when 
length and depth of successive pulses are. Ci) identical. 

Cii) different. i«ihile considering single pulse excitation, 
it is of interest to calculate response parameters i^ar the 
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Starting point of successive pulses. Sntih ^ inf Qrraati<Mi 
later on helps in deteriaining the response of vehicle dne to 
the excitation of series of ptalses. (b) Fourier analysis* 

Fourier analysis is an important tool for solving different 
type of excitations. Since these pulses occur periodically# 
equivalent Fourier series can be considered for the excitation. 
In the Fourier analysis response for each harmonic con^nent 
is calculated and then response of vehicle due to the 
excitation of Fourier series is obtained by superinposing 
response of each harmonic component. In Chapter 3# response 
parameters for harmonic excitation have been already calculated. 

In Chapter 5# length and depth of pulses which were 
deterministic in previous chapter are taken as randcmt. Such 
an analysis is realistic in calculating fatigue damage of the 
joint and vehicle components, A fatigue failure is the resiilt 
of cumulative damage that arises when response of structure 
to external excitation fluctuates. Uniform prch ability 
distributions have been assumed for random variables. First 
and second order statistics of response parameters for -fehe 
n^ pulse are detemdned and then Tof using superposition 
principle# second order statistics for series of n pulses are 
obtained# 

/ 

In Chapter 6# f i^t and second order statistics of 
the response parameters to continuous unevenness are calculated# 
The continuous unevenness has been treated as a hconogenous 
rar^om process with|Sero i^an# hater on continuous uneven%i .^§9 

5 

J- 



6 

is stiperinqposed upon discrete unevenness due to random low 
Joints, Second order statistics of the response parameters 
for total unevenness is than obtained by assuming the two 
processes are uncorrelated. 

Summary and csonclusions are described in Chapter 7, 

Following are scane of the in^rtant conclusions drawn 
from the analysis, 

(A) fi-esponse of vehicle to discrete unevenness shows 

that* 

(i) Peak values of maxlRum absolute displacement# 
acceleration and contact force occur for a particular value 
of dimensionless ratio rs<J.,l in our analysis). This critical 
value of r depends upon velocity# length of pulse# natiiral 
frecjuency and damping coefficient. Typical peek values for 
r m 1 , 1 # natural f retgpiency u)o *» 13,808 r&d/aec and |»ilse 
length w 1*5 bi are as follows. 

Peak value of absolute di^lacement is about 5,3 tjUmes 
depth of pulse# al>solute acceleration is 0,5 tis^s acceleration 
due to gravity and contact force is 1,45 times the total wei^t 
of the vehicle, ■ , ■ 

i3JL} Haximum. values of response parameters vary 
linearly with 2oiv (idjere is the depth angle » 4 titees ratio 
of depth to the length of pulse and v is the velocity) # keeping 
ttx and r as fixed. 
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fill) liffiitiiig hm. £0m& lOsovet 

veMoia it* contact to the gjpcwi* 1hi« tfelocitf iSf 

la^ur for t m i*l« 

fiir) fa) For a taloeitif leas than $6 Jsm/ke^ iateracticKiii 
effect ^previous jji^ee i» aig^nif leant* 

(b) for a velocity linf hmtwmm 5g to 10© 
effect of at leant 3 previoti# paisas alwuM he conaMerel* 

fe) For a velocity greater than i©0 Ica^ir* effect 
of atieast S pulaes ia neceseary in aoaaiatioii* 

M For a velocity of la0ir and a variation of 
pnlae length frora 1*5 laeters to 7*S aetere# range of 
( F^y/ fW ♦ «)«) la from 1*27 to 1*7«* i^le for <Fjjjj^jj/(M 4 »)g) 
ia ©*42 to ©*fS* a®d F^^ are oaae* iRin* contact 

force and fM e a©g ia tlie total veiglht of the vehicle* 

■ Probability of getting F|^^^Cii eii^g (alao hnovn m 
dymanic load factor {SitlP} greater than l*S© ia very hl©h 
f0«8)# vhile placability nf getting ©l*F greater than 1«7S ia 
only ©«1* 

fvi) haaunlng the design J3l*F .aa. 1*77 at a velocity of 
iO kat^iTf the prob^ility of failnre for one loa* length of 
tr.ick» ia aa high m It increases as ve (^crease tiM. 

value of design SltF.* 

fB) ^itesponse < 9 ^ vehicle to conti'naons unevenness shswi 

that 

(i) B*ai*s« vales acceleration contact §mmt 
.immmmmr linearly at high valocil^* Typical values at a 



8 


velocity of 150 Ica^/hr are “ 0*27 times acceleration due 
to gravity eaad *» 0.24 times the veight of vehicle. Results 
vary imieh depend upon chosen value of rou^ness constant, 

(C) Response of vehicle to total unevenness (discrete 
and continuous) shows that« 

1) Expected value of total response process is same as 
that of the expected value of response process due to discrete 
unevenness, 

ii) R.m.s, values of total response process increases 
with velocity* 


V 
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CHAPTER 2 

mathematical MOiDELE m TRACK AHL VEHICLE 

introduction 

A radllway vehicle is supported over wheels attached to 
the undercarriage system. The vehicle undercarriage system 
consists of springs and shock absorbers, supported on wheels# 
Railway wagons generally use Imainated springs* A sketch 
, ill, i.ig»igi>'' of the vehicle is shown in S'ig*2«i« 

A rail track as shown in Pig* 2* 2, consists of 
fishplated rail Joints and is supported by sleepers# The 
sleepers hold the rail in <x)rreet alignment and transmitiS the 
track load to the ballast* 

This chapter describes the construction of the 
idealized models for vehicle and traCk unevenness* Eguaticsis 
of motion are derived and transformed to a syst^ of first 
order uncoupled e<pation of motions# 


DESCRIPTION OP MODELS 


2«1 Vehicle Model 

The vehicle as described consists of several degrees 
of freedcMS and the analysis of su^ a complicated model is 
quite# difficult* Hence to sisqplify# it is assueed that there 
is only one degrees of freedom in vertical direction and 
idiole vdbicle mass (called sprung mass) is shunted on the 
linearized spring and dmiperii The mass of whe^el# aaele etc** 












10 


(called TiiD^ruiig mass) is concentrated on wheel itself. Fig, 2, 3 
shows the linear aKxael with single degree of freedoa. In the 
model M and m are sprung and unsprung masses respectively. K 
and c are linear spring and danping coefficients, Z(S) and 
Y(S) are absolute displacement of sprung and unsprung masses 
respectively. 

Vehicle Motjcm t 

Vehicle Is assumed to be moving with constant 
velocity v given by 
S a Vt 

where S and t are distance and time respectively. 

Parameter Investigated t 

Interest lies in the calculation of the following 
parameters, 

i) Maxifinm value of s^solute displacement Z(5} of M, 
li) Hsximm value of i^solute acceleration t{S) of H 
lii> Haxiimim and minimum value of contact force 
P(S), where F<S> is defined as 

F(S) « mU} ♦ ®YCS) + (M + m) g 
where <,) denotes differentiation with respect to time# g is 
the acceleration due to gravi%, 

tt,Z Models of Yradk Vmivenness 

Ihe unevenness of the track can be considered to be 


of two types# 


i) Eiscret* tineveimess to hogging of the rails at 
th« Joints, 

ii) continnons tmewnness due to the irregularity of the 
rail profile, 

2 * 2.1 nmmumm 

The shape of the discrete unevenrmss of the joints 
can be assuiaed to be of the following form (fig. 2. 4). 

TCS) m «hC 1 sin ^ * nL) ) when nh ^ S 4nL 

V -hid ->8in ^ (S -nil)) when nh S nl* + 
w 0 otherwise 
where 1 » length of pulse 

I* * length between two consecutive pulses 
n 9 nusber of pulse considered 
h m maximim depression or depth of the rails. 

When a vehicle passes over sueh a pulse# the response 
parameters are calculated by two methods. 

(i) Single pulse Excitationt 

In this analysis, response parameters are calculated 
for single pulse excitation only* Response for series of 

pulses is ^hen obtained by superimposing response for each 

I 

imtlse* ’l 

lii) yi^arier analysist 

Since the occurrences of the pulses are peric^ic. ; 
Fourier analysis can be applied here.* Fourier analysis.:''' is 
easy to apply for any pulse shape, egalFoleot Pov^ier 
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series can be ^iritben as 

y<s) - 2| * t , CO. * »>„ ‘in ^ 

» ^ I* Jt# , 

Ba»l 

where &, & ana are Fourier coefficients aiKl can be 
o n n 

obtained by following integralsi 
L 

% J if«s) as 

© 

2 ^ 

On • f j ifts) cos as 

1j 

^n ** b / sin as 

o 

Sino^ Fourier series is a ccxobination of sine and cosine 
terms# the analysis of sinusoidal excitation is necessary 
for calculating response of v^icle due to e&dti haxmonic 
component. 

The shape of the sinusoidal excitation can be assumed 
to be of the form 
yCs) as a sin 

where h is the wave length and a is the amplitude of 
undulations* 

In Chapter 3, sinusoidal excitation of the above form 
is considered and respcsise parameters are calculated* 

(^DtiiRtsus Paevenness i ^Flg*2.$) 

Continues unevenness due to irregularities of the 
rail profile has been assumed to be^ homogenous randta process 
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with zero mean and spectral d€nsity given byi 


§5nrm ) 


2 ; 

n. + Y 


j n 1 6 ao 


« 0 otherwise. 

Where C and Y are constants* £1 is spatial frequency# 
The constant C is given by 



where <r Is the R.m.s# value of unevenness. 


2*3 Bcfuation of Motion 

AS shown in Fig* 2. 3 the Inertia force due to mass M is 

# # e ♦ 

WZ . where Z is the absolute acceleration of mass M, The base 
is given a displacement Y{S)« The spring forge is due to the 
deformation of spring and is equal to K(Z • y). The daaping 

# m 

force is proportional to relative velocity (Z » Y> and is 

# • 

equal to the C{ Z Y). Suimning up the forces gives 

HZ *► C(z - ¥)+ K (Z * y> * 0 
# # # # 

or M Z + CZ + KZ » CY * 1 - KY 2.1 

where (*) denotes differentiation with respect to time 
dividing throughout by M we get 

writing ^ « 2 ^ U)o and ^ 

we get 

2irujo 4 t « a?u>o t •»- o)<?r 


2*2 
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Where ? is the dancing factor and Wo is the natxaral 
frequency. 

For constant -velocity run 



ss V 


hence equaticai 2*2 can be written in the space dosiain as 
+ 2Te uJ 6 vZ‘ + Wo^Z « 2 ?Ui 6 VY* + u)o^ 2.3 

where (*) denotes differentiation with respect to distance 9 . 
Equation 2.2 can be expressed as a system of first order 
equations (see Appendix 1 ) in the form 


dx. 

dt 


IT 

i^i 


I 2 ? U)o Y ^ 


2.4 


i B 1#2 

where o(^ w ..ij 

M. « 2 4>? (TgUlot «.) 

JIn> ^li* 

and Is 1^® oigen Yector of the dynamic matrix (Appendix I) . 
aljove equation cm also be written in space domain as 


dx 


<13 


% ^ t ujS Y^) 

% 


^ • S ** *1 * "W ^ ' dt ~ ) 


2*5 


1#2 


I 


For constant Velocity run 
dS ' 

a V* 


Substituting this walue equation 2.S becomes 
^ |2?u>oY*v t u)o Y I 

1 B 1 * 2 


2.5 
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2*4 Solution of EquatioEt of Motion 

Equation 2*6 is a consteffit coefficient first order 
linear differential equation. The hcanogenous solution is 
represented by 


*» e i i 

where are constants* 


Xi 2 


Substituting in hcaaogenous part of 2,6, we get 


VA. 


or 


ss 


V 


hence Xj 




" V 


^ih * h 
h stands for homogenous. 


i 1, 2 


2,7 


are constants to be determined from Initial conditions, 
Kon-hcMRogenous solution can be written in terms of 
convolution integral 


S 

*inh “ j *’i‘P> ^ 2.8 

where P^Cp) « ^ 

h,(S-p) - ^ -P > 

J, ^ 

also for constant velocity motion ^ v 

substituting these values# we get 

^ 2.8 

O ^ ^ i 

Absolute displacement 2(6) Is the stoKmation of hdraogenous 


and non«*h<»ogenoU8 part, 
2(S} a 


2 

X 

I«1 


Vih ^ 


a 

51 
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CK-S 
i V 


+ 


2 

2 

iasl 




S 


1 


& 



(S-p) 


{2?u>oT^ +U)fy|dp 


2.10 
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CHAPTER 3 

RESPOHSE OP VEHICLE TO SINUSOIDAL TRACK PROFILE 
INTRODUCTION 

Analysis of harmonic excitation to the vehicle is 
necessaxy as it will fonm a basis for Fourier analysis 
considered later in Chapter 4* Fourier series is the 
combination of sine and cosine terms. Hence for calculating 
the response of vehicle for Fourier series, response of vehicle 
due to each harmonic ooirponent is calculated and then 
superin^sed. 


3*1 Solution of Equation of Motion 


As shown in Fig, 2,6, sinusoidal tradk profile is 
represented by 
y « a sin 

stibstltuting the values of Y and Y* in equation 2,9 and 
changing variable from S to p* we get non-hCMogenous solution 


X 


inh 


J 


e -tr* (s —p) 

V 


^1 ( 


M, 


, 4 Y mova c os 2 tip 
K A A 


4 UJo 3 


2tr p 

A 


dp 


i 1, 2 


«iS, 


s a 


M 


® V f ( 4 Y u)o V Tt c os 2n p 

i o c A A 

^.u^3in ap 
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4)!s 

M, -J - ( M + Ulo ) 


J 

o 


.«jp 

V 


e ^ sin — *s + y j 


( 2 n p 

C } 

dp 


a. _■*! 4'?U)o«vn 

where y * pan 


*±nh * 


**t 


M 


V ( C 4tU)6Vn ) + u)o ^ ^ 

V ( > 


f e 

( e — 


-i" 




n4,j2llp 2i. cog { +Y H 

— sln^'^+Tj - -5P cos ( *^ )) 

( #4tn vu)a)^, ^ ( 

,v ( ^ 3 y z < 

t ( ^2 ) "3jy \ 


( 2 ns + Y ) - 2 n cos 5^ ^ ® + Y ^ - e ^ 
( T / ^ cos ^ ^ + n j V f 

«.nt. ^cosv}) 


Hcaraogenous solution as obtained in 2,7 

<H.S 

%« 


*ii> " 


3.2 


hence 


(i) M>solute Displaceiaenti 2(3) 


2CS) 


2 


z 

1^1 


*^1 *lh ■*' 


I *^l^ihh 

i««l 


<11* 
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Z H ^ 1 ® 

iasl 




1*1 
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( V 
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^ X 


( 21T (2 TT S _ ) - 

^ •» JbM cos ^ HhV j 
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V V ^ ) 


-. c( 


1 SinY - ^ COSY Y ) 


3,3 


ii) Absolute acceleration< 2 (S) 


4 Ih 


r CS) ^ J ocle% S g »i V } (l l j 4- u)o^ I 


Isl 


1=1 M, 


1 

<-«.) 2 
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ill) Contact Force * F(s) 
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^ q< 4S 2 a 4 > 4 *v 
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< 1=1 ^ ^ 


1 '^ 1 ® 

* e 4- 


M 


1 - *• 1 


2 C 

'^C 
C 


^ C 4t u)oyiT) ^ 4 ) 

( Xa ) 

( ) 


'^Ewr^ii 

nrn* f 


I 4 .Y| 4 ^p;os|< 2 ^ +y|| 

X ' ■ X 


«> 
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e 

i V 


(2 TT V) 
•ni( A 




a sin 


X 

X 

X 


V 


sinf - 


2^ cos X 


2 ns 


+ (M4ia)g 


) 

X ))> 
X )^> 
X n) 
) 


3.5 


2 l^snX'ts 

^^sults presen-ted in this section are for bogie 
CE- 28251, the data for xvhich is given in the table 3.1 
(Page 23) other data used in calculation are- 

i) Amplitude (a) » 0.1 

ii) wave length ( A ) =10 meters 

iii) Oamping ratio i%) =0.2 

iv) Range of velocity = 0 to 200 ki^br 

Following observations msy be made from the results 

(Fig. 3,1 to 3.7). 

1, Dimensionless ratio r(=2n v/AWo) is a very 
important parameter for the analysis, as it relates velocity 
to the wave length. Thus for any fixed velocity and natural 
frequency, increase in r iiiplies decrease in wave length and 
vice versa* Similarly for any fixed X and natural frequency 
increase in r implies, increase in velocity and vice versa. 


2. Ratio of the maximum value of absolute displacement 
to the amplitude (a) has been plotted in fig, 3.1 against r. 

It is seen that ratio increases very sharply at low 
velocities, peak value at f- = 1,0 and then decreases ^rery 
sharply upto r « 2,0. This critical value of r depends upon 
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dair^Jlng ratio ( value is about two times aR5>litu<3e* 

3* Ratio of the maxiimim valiie of absolute 
acceleration to the acceleration due to gravity has been 
plotted against r, in Fig,3«2 Ratio increases very sharply 
at low velocities# peaik value at r » i,G and then decreases 
slowly* Peak value of ratio is about 5. 25 times acceleration 
due to gravity* 

4* Ratio of the maximum value of contact force to 
the total weight of the vehicle has been plotted against r 
in fig 3*3 Ratio increases very sharply at low velocities 
peak at r =3 1*0 and then decreases slowly. Peak value is 
about 6 times the total wei#it of the v^icle* 

5* In all the fi.gures> it is seen that a particular 
value of r give all the peak values, This critical value of 
r depends upon natural frequency# velocity and danplng ratio 
and is important for the design of vehicle undercarriage 
system* 

6, Ratio of the minimum value of the contact force 
to the total wei^t of the vehicle has been plotted against 
r in figure 3*4* Ratio decreases very sharply upto r « 1,0 
and then increases slowly. It Is seen that at r » ,62# 
contact force is 'zero which implies loss of contact of the 
vehicle to the ground. Since ve'.lcle losses its contact# 
governing equation of motions will not be applic^le after 
that velocity,;. 
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Raticss of absolute aisplacement to tdhie wpJLltade, 
absolute acceleretioa to the acceleration <ane to fravitjf aisQ 
contact force to the total 'weight of the vehicle have been 
plotted against S/a in figures 3*5 to 3* t* Sis the distance 
along the wave length a § value of r is 0*2* It is seen for 
that for simsoidel exlctatims* ratios also varr slraisoidaly 
Cas expected) « Also rngximiai positive and maaelffluia negative 
values are same* *Xhus it is o^ncluded frcstt the anal^is that 
a value of dimensionless ratio r gives peak values of 
displaceiwjnt* acceleraticxi and contact force ‘Itesults 'very 
imuih depend* upon chosen value of afc®4itude* Also a 
velocity can be found above which vehicle losses its 
contact to the groundii 
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TABLE 3,1 

VEHICLE CR - 282S1 

load in tons 

1 Tare 

2 Gross load 

3 Carrying capacity 
^ Ajcle load 
5. Unsprung weight (m) 

6 Sprung weight on each wheel 

a) Under tare 

b) Under gross 

7, Weight on each axle 

a) Under tare » 1, 8 x 2** 3«6 

b) Under gross a 9.2I0C2** 18.42 

8* Total sprung weight (M> » <40.64 - 3,80) 

ts 36.84 

9, Average spring stiffness « 0.179 tons/iwa 

There are four springs on two axleSf 

each axle having two springs. All springs 
are assumed to be in parallel combination. 

Total spring stiffness « 4 x.l79 

m 0.716 tons/wR 

10. The nattiral frequency of the 
vehicle in determined as iSollowsi 

U)5 « 0.716 X 1000/ (36,84/9,81) 


=s 11.0 

a 40.64 

sst 29 . 64 

*# 40 * 64 12 sp 20 .32 
a 3.8 

* 1.80 
= 9.21 


m 190.66 

or o)©#* 13*808 rad/se<Si| 
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CliAPTER 4 

RESPCiaSE OF VEHIOiE TO PERIODICALLY OCOJRRING LOW JOINTS 

INTROEUCTIOK 

A rail track consists of fishplated rail Joint and 
is supported on sleeper and ballast. The sleepers hold the 
rail in correct alignment and transmitt''’ the track load to — 
the ballast over a wide area, ^ rail joint is the weakest 
part of track and a study of the effect of the Joints on 
vertical response and track forces is veary inportant from 
the point of view of design of the vehicle and track, A 
common feature of these Joints is the hogging of the rails with 
passage of heavily loaded vehicle, when hogging exceeds 
a certain limit, the Joints are called low Joints, The 
shape of the Joint or pulse (as frequently used in subsequent 
chapters) can be assumed to be of the following form (fig,4,l) 

y(S) SB • h (l+sin Si (S -L) when L « ^ ~ ® ^ ^ 

m - h(l - siniiCs -I*) when L s < L ^4 

£ 

ss 0 Otherwise, 

where ^ *s ^ and we have considered only one pulse. 

In this chapter first, response parameters for single 
pulse are calculated by using convolution integral approach. 

In section 4.12, response for series of pulses of different 
1 and h are obtained by considering any pulse and than 
summing up responses of each pulse,,,. In section a special 
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case of same 1 and h is considered. Since tJiese Joints 
occur periodically, equivalent Fourier series csffi be 
considered and response statistics are obtained in section 
4.21. Following parameters are calculated in eadh case. 

1) Mjsolute displacement Z(S) of sprung mass M, 

• ♦ 

ii) Absolute acceleration 2 (S) of sprung mass M, 
iii) Contact Force F (S) . 


4*1 Single Pul se Excitation 


Ecpaation of motion for single degree of freedom 
mo "el in space coordinates is 


v^Z* +2?uj^* 4- cOoSt «2tu)ovy*4- from equaticai 2-3 
Solution of equation of motloni 

Since y{S) is defined separately for three different 
region, we have to calculate response parameters separately 
for each region. Response for each region is calculated by 
using convolution integral approach. 


(a) First Region (Ii — ^^S^ L) 

Homogenous solutioni This can be written directly 
from equation 2.7 as 


•^ih 




where S. and o(- are same as defined in equation 2.7, 
X i 

Horwhoraogenous solutioni From equation 2,9 w® have 


^i»h 


J 


L -1 

f 


1 2 5u)ov ^ + wfy I <ap 
^ i 
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where Y as - h(l + sin si (p -.!»)) 


h^{S-p> 


hn. cos SI {p -L) 

OCj 

e -I <s-p) . 


substituting these values in above equation, we get 

iiS im ^ _ % 

^inh ® ^ V ^ ^(artUgVn cosO. 4- u)o 

J _ V Mj *' 

ft ^ ^ 

5 

( 1 +sln n {p -LJ^ dp. 

-h d> (1 . S ( 

” t ( “ *5"'' A (-U sin( n(S-L)+Y )- fi 

MV ( ( 

5 ‘ 

cos (0(s-.L) +Y)) + ~ ^u(S-L+ • ^ I -u4osY4“ilSinY p 4*1 


where 


C ^ -.1) 


m 4> 


r w 


fi 

V 

n V 
1U)0 


hence 


— Kz * ^ 

1 V 

2 

2 (XtOo-^ i » 1,2 


U)?(l +(2-cr)*>‘’ 


z(s) 


ZCS) m "E t 2 ‘^i^lnh equation 2,10 

±-i 



(. 


) 


^asin (n{S-L) 4'Y)4 -£IcosC + Y ) % 

( ) 

+ 'ucosy-asinv >’ 


( 


)) 

)) 


4.2 


Absolute Acceleration Z(S} 

2 2 oqS 2 V, . 2 ^ 2 1 

"S (S) « 2 2 ^ -u)o UCS-L+ |) 

^ss^. ^ M , # ^ 


i*l 
2 < 


i»l 


( u 


+ I ri^ (V -usiii(£i<s-L)+-?t-))-xi cos{a(s-L)+Y ) ! 

A ^ 

+ ^ U^e^^Cs-L-f ^ ) I acosY-^JlsinY ) j 


4.3 


Hi) Contact Force (F{S) 


F(S} sMZ(S) >f fiJSf (S) + (M + Bi}^ 




substituting values of Y(S) and 2 (s) we get 


( 2 
( 

C 


F(s) » H ( 2 ♦.j !lif + I !L^ „ r-^ 


i»l 


^ V ^ ( 




+ ^ J-asin( SI (S-.i.)+Y ) -noosCSl Cs-L) + Y )) 4^ 

{ ^ 
u?. ■*'|) ( UGOsY-^slBY^j+ m Wo r\ 

( » 

sin<SiCs-l»)) t (Mt ia)g 4*4 


(b) I S£ll S^I .,< I ) 

H<»K»genous solutlont 3?his can be written directly from 


eguation 3,7 as 


WIM 
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^on-hoEnogenous solutions From equation^we have 

ae , hj (S-p) / v j. niJ^V } 

M. 


inh 


j V + ujo^l 


dp 


Jj 


J . 


^(S-p) V + UJO^Y >dp 


where first integral has solved in equation 4,1 
For second integral 

Y as -h (1 - sin a{p -L) ) 

^ at hi! COS ( O <p - L> > 

0 <j 

and h^Cs-p) » e — vS- p) i a 1,2 


"in^ 


^ oC d> 

J - he — ;j~ (S-P) ( 2 ?U)oVSl COsQ-dP- h) 


ll 

2 


2 

+ ^sj(l +sin Q {p-L))dp - I - ( - h e (S,.p) 


V 


$•1 


2?WoVS2oosQ{p <^) + ojodl - sinQ(p - LJ) dp 4.5 

^4^ 1“^ |l - e -i^ij^lLcosY ^.usinY ) + | ® 


) 

) 


(aeosY -asiny) + ||l • - A ^ U.si« 


( 


I 


iais^h) +Yy^i).eos (a<s-i,>+vJ 1+ |e 

¥ t) 

*usinY 4 - Qcosr* ( 

* X > 


4.6 



hence 


i) Absolute displacement t 2(3) 
2 2 

Z(S) « ^ ^ 


i«l 

2 

2 


=^inh 


i.3si 




2 

4- 2 


^ 1«1 M^ir { 

U1 ( 


4>, 5 u)f? ^ I ' ^ ^ 


i i^l 


i i 


(1- e 2)) 

)) 


S ^ ^ T e i -nn ' ' 

4* A ( n.COSY 4- USlnY > + 2 (uCOSY - SI sinY ) 

( ) ^ ( ■ 


^U(S-Ij)j^ 2 ju8ln{ il(S-L)+Y ) +C1 

(n(&.L) 4-Y )} 4- T e !usinY4'SiCOSY H 

) A ( jj 


*-^ i l 

COS 


4,7 


m ♦ 

11) Absolute acceleration s 2{s) 


2 2 

Z* (S) « 


i*»l' 


2 ^ h V ( 2 

e — rr* +• ^ ■ „ i. .» , . JL.i,. . ..»i. . /-'uib© 

t»l ( 


u{S-b) I ii2 (« asin(a(S-L) 4-V )-.Q:cosCQ(S-I,) 
A I 


Y )) ) + T {^sinY+QooaY) ) 

h ^ < )) 


4»8 


iii) Contact Fotcet F<S) 

P(s) » M Z{S) 4-t«i*(s5 + (M 4- lB)g 

# ♦ • » # 

substituting values of Zi(S) and Y{S), w® get 


X 2 * 2 . «1® 

F(S) a* M I 2* ^4© 


t ^ 


h ^ V ( 


1-1 M. 


(-ClUJo 

C 


B ^ |••u.siIl { S2{S«»I<) +Y ) wSic^s { S2 (S*.!*) 4*^ ) 

) ^ ^ ^)X 
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(c) Third Regions {S^l»+ 


1 

2 


) 


In thisregion, there is no ground induced excitation 
since Y{S) is zero, ^ence this part contains only transient 
vibration which dies out as distance increases. The 
response for this part can be written in the form. 


2(S) ■ g(S-S^) + Z g^h(S - S^) 


hence# 


i) Absolute displacen^ti 2{S) 


2(S) « Zg^ g(S Sq) + h{S - sj 


So 


4,10 


where L + — 

o 2 


So 


2 * . ^ Jr mi 

» Z(S ) « y ^ ^ ( 

0 ^ mmmmm 




, ) T 

^ ^ ^ 1 - e — I j 


( ^ T U1 ^ ^ U! 

( Q. 000^4* usinr v + T ® (uoosY-SlsinY ) + — 


ucosY-aslny J - f 


X T U1 t 
* A ®*2 


* 2 *1 ^ > 

i =» 1 ' 

2 ^ h 

a (s^) . 2 '»>i \ e <1. + i )+^2^ ^ 


^So “ 


iasl 


^oTX 

^ j - u sinY-42cos y 2** U>o® A 
t )) 

C U sinY-tOcGSY 

< » 


2 

ul 
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g(S -S^) at e V %■ jCos g-CS-S^) + -^£-->!£y^-- 
Sin'SS (s _s^) I 
h(S-S^) - ° 

B u)o<i 


4.13 

4 . 14 

4. 15 


il) Absolute accelera'tlon* Z (S) . 

“ ^So’S + ^so 

• • ^ (S*»S d)<*? 

where §{S»s ) « e v ® ( sin (S—s ) 

O ^ w o 

X Wo^- ^«3^)+2?U)oWdX+ cos 

Cs-s^) |,^».I. (-2ru)oW^) +t;^U)o^- 


( 

h (s**S ) » 
o 


- {s-^ ) 2 2 .2, , ^d 

J (X Wo - ^<5 ) sin — 


is^sj - 2XU)oWd Cos <S-S jJ 
O VO) 


4.16 


ill) Contact gprce ? F{s) 

F(s) » M Z (S) 4 - (M+nj)g 

« M(Zg^g* {S«s^) 4Zg^h*(S-S^»4- (M+i#g 4.17 

4.1*2. Excitation di^e to Se ries of Pulses 

In the previous analysis of single pulse excitation, 
we have calculated response paraiaeters by considertng 
single pulse excitation only, When series of pulses of 
the some shape best of different size (different 1 and h) 
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occur periodically* then response is calculated by 
considering any pulse and then suRaaing up response 
du€ to previous pulses <By using superposition principle) , 

For any n^ pulse, shape can be described as {Fig, 4. 2) 
y^<S) *-hjj(l 4 -sin fin{S -nL)) when nL- -| ^ S £ nL 

K 

sat • h^i 1 - sin51n{S-.nI»)) when nL 4 S^nL + 
a* Q otherwise, 
where Sin *, 
n « 1*2 — 

In Fig* 4*3* and are length of pulse and 

®th respectively, ^n depths of 

**th P^^®® respectively. 


Response parameters for any n^ pulse can be 
written directly from the previous single pulse excitation 
analysis. Hence for three region, response pararater will 
be as follows s 


(a) First Regiom nL 


4 s ^ nit 


i) Absolute displaceB»ot it Z (S) 

X ^ ^ 

Z„ , (S> m E <S-.ni.4- Jbl )+ 2 \ I 

X r V 2 ^ „ i 

t C 


(1- I usln(Qii(s-id») 4- Yn) 

II ' 

4- fiin cost ClnCs-ia,>4- Yn) |4- « u.(s.-ai:«4' 'tst t 

^ K 2 

I ucosYa -Sin sia 4»18 
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where Is the absolute displacement due to 


pulse only* 


T 

n 

m 

X C2? U)o tfij.) ^ 4 COo X 


m 

u)q X 1 4(2"^ ^ ^ 

^n 


nv/ljxio 

n 




■Y„ 

n 


tan*^ 2 '? 

u 

set 

\/y 


11) Aljsalute , ^cpeleratlon. V (S) 

^ 2 oi 1 ^ 

Z -CS) « y a oc e (S -nL-f— §3 + 
n, 1 ^ *^1 i i V i* 1 

i*l 


V ( 2 u(s-nL + •*“n) T 2... X 


n ^1 

1C 


( 


u oDo* ® 


2 + Ji X 

n 


h 

-usln ( £ln(S- nL) + Yn) -Sin cos (^«(S-nL) | 

1 ^ 
+ ^ ^'J-(S»nL + *-|3 Xacos Yn -On sin Yn % ) 

n 


%) 


where 2^ (S) is the acceleration due to pulse only. 


iii) Contact Forget 


^ + Cm 4^ hO 


2 2 c<4 2 

X *=<1 4> .e -§.{S-nL -1-3 ) ^ 


2 ^ -3 < .u ujo^e j-5 + 


C 

» mC 
C i 

2 
1 
i « 1 


.19 
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>'-usin {£ln(s -nl.) + Tfn) oos (On (s -nL) +Yn; 


X . — n ^2 a(s^j^ .jH )^ucos’Sin - Sin alntn i + -^ U, 

2 X * ^n Z 


X^ 2 A 


X^!n 1.2 


n 


u<s-ia 4 + n ) 


2 *5 ucx»s Yii -S5n sin j + m u>o‘ 


h^sinCln (S-nL) + (M g 

H H ^ 


4.20 


n 


(b) Second Reqion t nL£: S £: nli + — ^ 


(1) Absolute Displaceinent t Z -Cs). 

** § m 

3 ^4 2 ^2 
z„ ,<s) m Z e. <t>. e --*<s-ni4) + S ^1 

^ Issl M^V 

V 2 uX SC ) a? 

5 ^^0 (1 - ® ^ ^ < Q^s vn +usin Tn > 4> ~ 

J U n < ) n 

(ucos Ya -fij^slnYa ) + (1 -e ^ 

z u «jj 


Y ^ 

c vJL sln (QnCS-nl.) +'<‘n) +Qn cos { Sln(S-nI.) -► Tn) + 

* X ” 

^U(s-ifli> Yn +Q.n cos Yn I | 4.21 

where Z^ ^(s) is the response due to the n^ pulse only. 


♦ • 


11) Absolute acceleratlcm f 2jj^^2<S>, 

. . 2 2 ^4 ^ W 

^(S) s* ^ Pj^e (S-hL) + ^ ®||' V 

i.i * 


a, 2 


1»1 


( 2 U(S-nb) ^n o 2 ^ 

UUJo® • — 


2 ' 

n ( usln (Q.a{s - nl») 4‘Yn)-^n 


€SO» 


{ 


a 


^n{S- m.) * Ya ):)+ ^ e *“ 3S UsinCln -f^^n cos YU) 


4.22 



3S 


I4iere 

• m 

acceleration due to the pulse only. 


iii) Contact Force t «(S) 

iw«iiT-.Tiriirrniiiiiiirrtirniriini:iiiii-iti iiiiii;riinnni |■,rnlI■lB^llr ^ 




# • 


2 

i« 1 


*^4 

n 

M. 


s M< 2 "^(S-nL) + S 

( i* 1 i ^ ^ ^ i« 

f ■» uuio& ^ f — usin CSlnCs— nil) 4- Tfn) 

{ n ? 

-iin cos(Q^n<s-.ni:#) tVn)j + 

0 

^ UslnYn +0.n cosYn|; j- ra uiot^ slnSln(S-.nL) +(M4*a)g (4,23) 

( ^1) 

Wheie contact force due to the n^ pulse 

only. 


l^hird Region t nL 

(i) Absolute dis-ol..sgement t 

^ (S) « Zjj.^g(s-sn) t ig^h{S-Sn) 
i 

Where « nL -i- n 
n y- 


4, 24 


* 2(S ) 
Sn n' 


U .1 


I ^n X 2C0o^ ( , . _ 

1 i 1 Itv X U ( ^ " ® 2 

X ' 

) % 


n 


T ( ) T » 1 1 

■*■ - l^n cos Yn+usln Ynj-j- (ucosYn '*»Cln sin Yn ) (e 1) 

a f ) n 4l 


^ ^ ^ i 

e (usin Yn -i-ftn cos Ynj|+ ^ *^^4 Pjt® 

*n ~ ( )i i « 2. 2 


Ttr, « z(s ) 

Sn 'a 


I «i f. .c-a i|,^ 

i « 1 ^ a V 4^ 


2 

4,25 


2 

Z 

i. ssl 
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2 

^ X X 2 ul« 

{(«- "■ ■' -£% X^usjimifn - Qn cos '^n X - u)oe •r;— 


"i ( n X X 

* li3._ f O 

— e— ^ ^usinrn+Qn cos 


4*26 


n 


g is *» 


sin (S . 3^)1 
or-/l gi(S-.s^) s ® 


I cos ^ Cs-Sn> + 


„ (s«.s ) 

V n 


4# 27 


(^os (S.3J + 'S 

c ^ 


) 


and h (S-s^) =* 

U 


. _ sin^ (S-Sn) ) 

(1 ^X^)h > 

<S.S ) 

V li' 


eJ^ (S-S ) 

^ ^ sin 


coa 


4.28 


tiJd » 0)o <1 - ’5^)% 


Where 3^35 is the response due to pulse only. 


w w 

(ii) Absolute Aceeleration t Z^ 3(S) 
'\,3(s) ‘^sn9 ‘=-=n> ^ 


4.29 


vdiere ^ 

... _ . -2Lii>P(S.S ) if^d /c._p ) X^X__„ 

gMs^S > « e V n ^sin — xTi- 7 ^ 


n' - > » 1^ (1- 

(?*U)0 -wa*) + 2iU)6'Sd l^cos ^ (S-S^) 

{ - 2tU)o>Od) +lJ*COo* 4.30 

Wa 


• * s'* (S-S ) , 

ll CS^S^) «« ,■■■.,, c:,!--. ,■-■--- ^ sin y 


(S-.SJX%^U)0^ 


J^d^ j[- 2?u)oWd <30« 4,31 

# « 

where Z^ ^is the acceleration due to pulse only. 
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(iii) Ckantact forces 

V3<S) = “ ■™’® 


MCZg^g(s-Sjj) + Zg^h (s«Sjj)) + {m •«ft)g 


4.32 


Here# above# we have calculated response parameters due 
dxie to the pulse only. Response parameters for series of 
pulses are obtained by summing response parameters due to 
each pulse (using superposition principle) , 


Lute d-lj 


Z(s) 


M JL 

zm « I "li# 3 ^^^ "*■ ^14+1# CN-fl)L - 


+ 1 ) h 

II 2 ,. 

-Z ,(S) -KS„ . -(S); when S (51+1) I, + 

K +1 1 


a ^ z;„ •aiej\ When S^(N+l)L + - -"p -g 

n « I 2 

where N is the number of pulses passed* 


4,33 


2) Absolute acceleratlom Z(S) 

** 1 


S 6 (N+DL 
N 

^ ^ ^a,3(s)’*' ^+1,2(3)? (S+1)L^S^ (N+DI. + -y 

1 

5l’^n.3(S) When S J(H+1)I. +-a|l 


4^34 
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3) Contact Force s PCs) 

PCS) a m'zCs) + CM + m)g + m Y* (S) 


hence# „ 

C ^ 

FCS) aUC^ z 
Cn=:3 

when CN 


) * • 

n,3CS)'‘‘ ^+i,iCS) ) tCM«R)g -mr^is) 


( ^ .. ) 

° ” (nrf ^'■>.3(S)+"^!+i.a(S) I * J«‘y„<S) 

1 

men S < (H+DL + 


(K+ u. ) 

Vacs) 

' Ilsrl ^ 

X 

When s $ CN+1) I* ■<- 


4.35 


4 . 133 . Special Case 


A special case can be considered when all the 

pulses are identical (same 1 and h) • In that case pulse 
shape can be described as (Pig. 4.4). 


Yjj Cs) *s • h C l+sin SI (s«>n!#))t when nL-» ^ ^ S 6 nli 

1 

- h(l- sinil(S-ia»))# when nL£: S4 ril. + - 
* 2 

m Q Otherwise# 

where SI «* ^ and n » 1# 2 — 

Response paraaneters for different regicn can be written 
directly from the previous seirtion 4.12. 
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(a) First Region $ ■T'L - -:S.S £ tdi 

2 

(i) Absolute Displacenient i 


'n,l(S) 


M **# * 

2 ^ *^i i ^ ^^u.(s».nL t 2 


i as 1 ^ 


^ A ( ^ Sin ( O.Cs««bL) •^ 'y)*f'Sl.GOs C ^ CS-«iiX») ^ 


- <6 

g'^{S-i3ii4- - ) yUcosY -Qsinv y I + - *^i 

<£ A A J X S3 X 

X X) 
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ii) ■^solute acceleratlou t 


• ii 

^n,l<s) 


2 /N ^4 i Ik 4* 

1 h tj; e <S -I* 4 - 2 ' **■ 


isx 1 
2 ( 


2 ^ ^4 2 u{s -m, 4 X ) f SI 

5r ± VC -uoio ® 2 + A 

* ^ < M> ( 
i SB i X ' 


{ I 

(- US1»( QCS-ni.) +Y> -ncos ( O (S-PL) 4 y) ) 


< 

+ S cosY-asia't'ij 


) 


) 
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iil) Contact Forces *’n#l(S) 

(S) t a'^niS) 4 <li4«i)g 


“ “ X 2 • 

X 4 SB i i 




(S«,ni44 2^ ^ 0 


h 4> 


i SS 1 
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( -uuooe 2 f {Sl<S-.nL) + Y) 

-acos (Q(s-.nL) +Y’)] j- I u^e uCs-nL-f | ) : 

D X 2 2 . 

XucosY- QsinY X ) X -t- ®tOor^ sin Q.{S-.nL) +(M«)g 4*38 

) X 

(b) Second Reaion A nL -l^S^nL 

1 

i) i^solute DlsPlacementt 


^n,2‘=> 


( uv^ 


“ "<1 ‘ ** 
2 V V I sJbt 

1 « 1 iol 


Y q* T Ul 

1 - e 2 + - {Qc^sY +u.sinY) 4- e 


t, 

J a x^- 

jUCOSY— SislnY X"*’~^^ (1 » e U(s-nL)) - ^ (^usin 

<a(S-nL) +'f>4.acos (Q(s-nL)+r))+ I 

) 

(usinY +Q.COSY )) 
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• « 

li) Absolute Acceleration! ^n, 2^®^ 


^ n,2<^> 


2 2 


2 't’i “i P.e -i (s-nL) + 2' 
iwl * 


2 

h H 


H. 

i=l i 


f -...A 2 i-tisin CaCS-nL) +Y ) - 

i -Ulbo ® A ( 

' ) 

a cos C a(S-nL)+Y)^4- f e^^-’^^JusinY+flcosY J) 4*40 


lii) Contact Force* 

^ ' 2 2 ^ 

» «( S <t>. “2 P.e <s-nl.) + 5 

( iil * ^ i= 


b 


C 

V ^ 

^ c 
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X-UWo e -.|Q.^(-usin +t)--Qcos { 

X 

Sl(S.nL) +Y) + |u e u(S«nL) X |0 

X 


in cOo^r^ h sinS^Cs-m,) + (M 4ia)g 
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(C) Third Region 8 S ^ nL 4. 


1) Absolute displacementi Z -(S) 

TX^ Jf 


V3(E) = »><® - S”’ 
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febere * nli + * 
n 2 

2 2 ^ 

^an - " Z *1 2’+ Z 

i. s 1 iasl 


'i ') t m 5|il 

JiM J 1 . e ^ A ^ 

3 u \ 2 A 


( a cos r-SL sin'Y ) + 


U)(? 


Xj. X JucosY-^lslnY I - ~ ^ 


( a sinYt Si COSY ) X 

X 

2 

^sn - ^(sn) 


2 '*’i “'l %® ^ + 2 


2 

i a 1 
2 ^ 2 


2 

%...±JL 
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< |-usinY-S2cosY |-a3? e^+¥e^ 

) 

{ uaiiiY +Sl<x>sY )) 


4,44 
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•I ^ U)d 

g(s-.s ) s e V <s -S ) (cos {s-s„) + 

^ H B ^ If , H 

g , sin (s -S„)) 

irt-pr^ “ ) 


h{s-s J s 

B 


(S-S ) 
•V’ n 


sin ^ te -s„) 
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ii) Absolute acceleration* « (S) 

^ + 2g^ h(s-s^) 

•• 1S> c 

where g (S-S^) * & v (S-.S ) ( 

* n n ^ 

+ 2?uJo>Ja)Csin (S-S^) 

+ X < - 2?u)j>aa) + -uJd^ Xcos 
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{s-s„) 


h' ) 
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r (S- S^) « 


^0 - 2>;u3^<3 oos'^ < S- s^) ) 


** (S^8 ^ ^ 

e V n' ^ fUd 

< sin -I (S-S„) 

|2s _ i C_ C \ ^ 
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(iii) Contact Force* F , (S) 

Miff ^ 

V3(S) ‘ " ^3 + <” -*>9 

e ♦ 

(* M(2g^g* (S-SjjJ + (M4ai)g 4*50 

t 

These a«e the response paraiaeters in differenfws sections Sue 
to Bii. pulse only, ^sponse parameters due to series 
of pulses are i^tained in similar to previous section. 
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1) -Absolute dlsplacement. > Z(S) 

^+1.1(S>' - 2 - 

H 

* I” Vacs) <N+l)L^ S^.(N+l)L+| 

Hal 

N+1 


2 ^n,3(S) v:hen S $: <K+i)L + ^ 

Ussl 
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2) Absolute Accelerations 2(s) 

z’(S)« f when (N+l)L- I^S ^ (N+1)L 


Ussl 


N 


^ (N+1)L£ S5(1?+1)L+ I 


n»l 
il+1 

^ *Zjj 3 ( 3 ) S^(H+1)L+| 

nal ' 


3) Contact forcS t F{S) 


F(S) « M 2(S) + <M4iR)g -im y^(S) 


( ^ 


> 


M 


( 5 2„ o<S) + 2 (S) )f (M4Bl)g 4mY„{S>,(H+l)L-|^ 

( Hal n+1 } n 2 


{ ^ 

% # < 


) <N+l)b 

a M ( ^ 2 ' (S) +’L^. «<S) ) + (M 4®)g 

(na 1 ) 


■HByjj(S) when CN+l)b^S^ (K^l)!- + | 

K+1 .• , 

» M( 2 : ,(S) ) + (M«)g J (N+1)L + • 

Hal 2 
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4.14 • Results 

The data used for the calculation is for sane bogie 
CR-.28251 given table 3.1. In addition to that data, other 
data used is given in T^le 4.1. (Page 55). 

Following observations may be made froai the results. 

1. The dimensionless ratio r (= itv/ltOo ) is an 
important parameter for the analysis, as it relates velocity, 
natural frequency and length of the pulse. For any fixed 
value of V and tOo decrease in r implies Increase in pulse 
length and vice versa. 

2. Response histories ha've been plotted as a function 

of distance alongthe trac3c in figs, 4.6 to 4,8 for r w 1.1. 
'R-»tt©5 of absolute displacement to the depth of pulse (h) , ^ 

absolute acceleration to acceleration due to gravity and 
contact force to the total weight of the vehicle have been 
plotted against S/1, Fixed parameters are (J0owl3,808 
ra<Vsec. 1 s* i.Sm, h as Siam, and L « 13m, 

There are two regions. First region is over the 
Joint (forced vibration) , and second region is beyond the 
Joint until next pulse starts (Free vibration) • Particular 
interest lies at the starting point of the next Joint, 

Maximum value of absolute displacement occurs near 
the end of joint \^ile the maximum values of acceleration 
and contact force occur at the centre of pulse. After the 
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pulse^F all parameters decay ej^nentially. For a Joint spacing 
of 13 meters, response parameters* at the starting point of 
the next pulse Is several order of magnitude less than the 
maxiimim value of response parameters due to pr^ce^dlng pulse* 
The interaction effect (response values at the starting point 
of successive pulses) is quite small* 

3* Ratios of maximum value of absolute displacement 
to the depth h of the pulse, maximum value of accelerations 
to acceleration due to gravity (g) , and maxiimm and minimum 
value of the contact force to the total weight of the vehicle 
have been plotted against r in figures 4*9 to 4*12* Fixed 
parameters are a)o« 13,808 rad/sec* h » Smm, Range of JP is 
0 to 4*4 (corresponding velocity range is 0 to 110 Icnylir 
for 1*5 m pulse length), 

4, Maximum absolute displacement increase^^ upto 
r at 1,1, then reduces raonotonlcally as r is increased, 
Amplification at r * 1,1 is of the order of 5,4, 

5, Maximum absolute acceleration increases upto rwl.l, 
re<&ices sharply upto r » 2,0, then Increases slowly as r is 
increased, Pei^ level at r » 1,1 is of tl» order of 0,5lg, 

6, Maximum contact force increases upto r » 1,1 reduces 
far r = l,l upto r 1,8 and thereafter increases almost 
lineariiy, as r is increased. Peak value is 1*46 times the 
total weight of the vehicle. 
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7. MiniitMHi contact force gives peak at r « 1*1, 
increases for r > 1*1 upto 2*0 and thereafter decreases 
almost linearly as r is increased* Peak value is 0*68 times 
the total weight of the vehicle* 

8* Peak values of maximum absolute displacement, 
acceleration and contact force occurs for r » 1*1 (this 
corresponds to velocity of 22^5 k%^r for l*5m pulse length) . 

9. The angle 2oc is also an important 

parameter for the analysis, as it describes the shape of 
the pulse* 

10* Maximum value of absolute displacement has been 
plotted against 2'o< v for different values of r in fig, 4* 14* 
Ratios of maximum value of acceleration to the g, and 
maximum and mininxim value of the contact force to the 
total weight of the vehicle have been plotted against 2«.v 
for different values of r in figs* 4.14 to 4.16 Range of 
2 oc V is 0 to 2.0 (corresponding velocity range is 0 to 
200 knv'hr). Other fixed parameters are 2o< ** 0.01334 rad. 

0)0* 13.808 rad/sec. 

11. for a particular value of r and2o<, increase in veloci 
higher 1 as r is fixed, but since 2oc is also fixed, 
incare ase in 1 increases h also. 

12* All the re^onses parameters vary linearly with 
2«<v* Maximum values oecair for r *1*1 Range of interest 
of 20 CV is generally 6.4 0.8* For r « 1,1, the maximum 

response values s in this ransp are given in the follo%dng 
table*. 
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“TABLE 4-2 


Response values 

2 0CVr! 0.4 

2c< vst 0.8 

Z 

max 

30mffi 

58 mm 


0.55 

1.10 

V 

1.52 

2.20 


0.65 

0.29 


13* F is called dynamic load factor (DltP) 

and it is seen from above that its value can be as high as 

2 . 2 . 

14. In fig. 4.16« for r * l.l, minimum value of contact 
force goes to zero at a velocity of about 80 JciVTir. This 
velocity is important as vehicle losses its contact to the 
groiind. For r » 1, 2» limiting velocity is 105km/hr, 

15. Figs. 4.17 to 4.19 have been plotted to determine 
the absolute displacement, acceleration and contact force 
of the vehicle due to excitation of one pulse, near the 
starting points of successive pulses. This interaction 
effect is important vdiile determining response parameters 
due to series of pulses. 

• # 

2^.. 2^^^^ and F^^^^ are maximum values of 

omax omax omax 

displacement# acceleration and contact force on the pulse 
n «0, the excitation effect of which is to be considered 
near the starting point of successive pulses^ 
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to the masclimsm 

value of displacement, acceleration and contact force near 
the starting point of successive n pulses (n »! to 5) « 

Displacement ratio (2 ^ ) acceleration ratio 

OIuclK.^ 

, ) and contact force ratio 

owlbx 


( (^t" ^ > 1 Plotted In 

Figs* 4.17, 4.18 and 4.19 against n for different values 
of r on log scale. Fixed parameters are tOo* 13*808rad/sec. 
L *= 13m. 1 xs l.Sffi. The reason of defining contact force 


ratio in the above form is that tlie value of 
<M-«a)g is nearly equal to 1, By substr acting 1 from 


and dividing it by (M+m) g) •!), it 

beccMnes easy and appropriate to draw it on log scale. A 


plot between the ratio of to F__ against n does 

*' nitigpc Uwoji. 

not show toy significant variation on log scale. 


16. All the values decrease exponentially against n. 
Also it is seen that as r Increases (velocity increases), 
all the ratios increase as clear from the following table. 
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17, Displacement# acceleration and contact force ratios 

have been plotted against velocity for different values of 
« =1.5 in figs* 4,20 to 4,22. Fixed parameters are k » 1.1# 
13.808 rad/^sec. I« = 13m. and 2o( s^. 01334 rad. All the 
ratios increase with velocity. Also at low velicity# even 
at the first pulse (n»sl) # ratios are very small, A clear 
picture can foe obtained frcan the numerical values given 
below. 


Velocity 

(ku^hr) 

n 

Di spl acemen t 
ratio 

Acceleration 

ratio 

Contact 

force 

ratio 

30 

1 

0.0208 

0.0233 

0.0228 

140 

1 

0.6280 

0.585 

0.575 

30 

6 

4.1x10"’^° 

IxlO 

0.9x10“ 

140 

5 

0.01127 

0.01246 

0.0120 


It is seen that at low velocity ( 4 SOknv^r) # 
interaction effect is alriKSst negligible at the 5th pulse, 

18. Displacement and acceleration of the vehicle due 

to the excitation of 5 successive pulses are the summations 

of displacement and acceleration due to individual pulse, 

evaluated at that point. Contact force due to excitation of 

5 pwlses is obtained by following relation, 

5 

S’*™ • *4 ^ + (M-m)g 

Here also it is appropriate to. subs tract from 

F ^_and divide it by {M4rai)g, 
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5 

^imax * ^ (Max values of 2(s) near the 

zi si 

starting point of 6^ pulse). 


^omax ** value of displjacement on the pulse n s»o. 

5 ' ' 

# ♦ • « 

^niaax ** ^ vahe of Zjj^ 2 near the 


staring point of 6^j^pulse>, 


2^™^^ * value of acceleration on the pulse n *0. 

oniax 

value of contact force on the pulse n asO* 

oiiax 


# • # • 

z to z z to z 

Umax omaxt tmax om^ 


Ratios of 

and ( “D^^^ve been 

plotted against a for three velocities in figs. 4.23 to 
4*25. Fixed parameters are r =1.1# uOq* 13,808 rad/sec, 
L 3 13 ir. 

It is seen that- 

i) For a velocity less than SO kii^/hr, interaction 
effect of only previews iwlse is significant. 


ii) For a velocity lies between 50 to 100 knv'hr* 

interaction effect of at least 3 pulses should 
be considered, and 

ill) For a velocity greater than lOOloVhr, effect of 
at least 5 pulses is necessary in summatiam. 


19. Ratios of maximum and minimum value of contact 
force to the total weic^t of the vehicle have been plotted 
against length of pulse for three velocities in Pigs, 4, 26 


and 4.27, tength of pulse varies from IfSm to i7*.5«i.“-r'- 






,r?'0537i 

6^- Mfl). V. 
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Corresponc31ng variation for h is Smm to 25inin* These ratios 
have been plotted to determine distribution functions of 
when probability density function of 1 is 
uniform. It is seen that the increase velocity increases 
pe^ value of max. contact force. While reverse effect is 
seen in the case of min contact force, 

20, Distribution functions of and 

(M+m) g have been plotted in figs. 4,28 and 4,29# for 

a velocity of Range of ^ 

1,78 and i‘^j| ^ (M4m)g is 0.42 to 0.95, It is seen that the 
probability of getting greater than 0,50 is 0,80 

while probability of getting P^^ ^CM%tt)g greater then 0.9 
is 0,05 only, 

21, Fatigue calcailatlons i Fatigue failure is the 
result of cumulati^re damage that arises when response of 
structure to external excitation fluctuates. As an application 
of the above analysis to the fatigue failure of the Joint, 

let us Consider one km, length of the track, the number of 
joints or load repetitions will be 

- 1000/13 - 77 

Since joint spacing is 13 meter 

Let us assume the design dynamic load factor (DLF) 
as as 1,77 for each type of unevenness. * ' 



(D1»F) can be e:Kpressed In terras of error function 


For contlnuons unevenness case# the mean value of 
contact force is the weight of the vehicle. Fa (M-i-in) g. 
Also r.ra.s. value of DLF at a velocity of 60 kr^/hr w Fa 
0.1375 {Mtra) g. 

fhere fojpe - 

P(DLF 1,77) a Probability of DLF or / (M+ra) g 

less than or equal to 1.77, 


1.77 - ( _ jy (Mtm) g) 
( y/(M+ro) g ) 


aO.S + erf 


1.77 - 1 
0,1375 


» 1 « 00 

hence probability of DLF 1,77 (probability of failure) 
is sero. 


(iii) Total Unevenness * 

Superimposing the two types of unevenness# the 
probability of design load factor greater than 1,77 
a probability of DLF 1*77 for descrete unevenness + 
probability of DLF 1,77 for continuous unevenness. 
m 0,96 + 0,0 
a 0,96 

Hence probability of failure for design DLF » 1,77 is 0,96 
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^ 1 


' TT ) 

C-h(l -sin -|* <S«L))) cos 

( ‘‘■V 


M f 

I 4 3 


{ S + cos sin (S-L) ) 


£ 14 . i 


2 f I 

j QQg J 2 Piiin.S — cos sin (S— L))dS 


1 <‘ 




jj, («in 2nna 

Mi mSmmm m 

^ ( 2 nTT>L 


_ 1 <cos n(S^ + ^ ) cos Tt(%ii -^) 




n t i _^ ) 

1 L 

) ll 

J ) 

) ) 

} ) 

I 1 


{ a>s tt ( s-l 4. ^-E > 

+ Js < ^ I. 

‘ IT< j + ^ ) 


COS n ( ) ) ) 

TI(| - ^ ) ) > 


I, + i 

X 2 


On solving# it gives 

. . ib X ^^n - "in °nu ][ 

^ ^ I (1 • 4n2^2JJ^ ^ 

I* 

^ ^ 2^ [ Y(B} Sin dS 

11 I* J " 


4*59 
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2J01I as 


^ j (14-sin -2^ (S-i,)) sin 

,2 

Ii4> ^ 

j (1 . sin (s-D) sin ds 


( - cos 2n TTs 
- { L 


i> ( 2mT/L, 


/S^I» 2nS V 

2 X 


X T1( 


2n 
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sin n( 


S-L 


2ns 


a. i?sa:.a: ) X ) 

r !■ 
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nc ^ + 


X) ^ (-<=®s ^iiis jg X»inTT(^ .2ns) 
2n ) X ) • L ( — i I — 1 ^ 

mmmm / % . >TT’ Jt « 


) , ( anTTA 

L-l| 


¥ n ( 1 „ 2n ) 

1 IT 

L-f ^ 


On solving it gives 


slr.n< 2^ + ^) * > ’ 


b^» 0 
n 


TT(|+ 


* i 

X |L 


Therefore Y(S) 


a «>o 

o 'c- 

-r- 4- Z a_cos — ^ 

2 n*l ^ 


2nTTS 


4.60 

4.61 


where 




2h U (1 - ^ ) 


^ ( gin bng* 2 U.n) 

and ^ ( 1 - 4ii^ )ii 


4.2.1.(b) Calculation of displacement, acceleration and 
Contact force! 


is 


Equation of motion for single degree of freedom model 


'z +2Xm^ t « 2 ^U)oY 4- uijy 


or in space coordinates 
u)o 


Z- + » 2? 
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Homogenous Solution 

Homogenous solution of the left hand side is 
ot.S 

^ih ** ^1® "v equation 2.7, 

(ii) Mon*homoqeno\is solution 


From equation 2.9 v® have 
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h^(S -p). 
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00 
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\ > 


4< 
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where 


I. « « ( 1 • e — j 


‘1 “ j 


4*63 


( . «^»<2a.H.§,. -V. ) 


<“ (2nTT)2 5j^( 

{ v)2 < L )2. ) 


^ e — j 


+ sin Y. ) 


4* 64 


X( V )2 nr-)2 J ' 


-V 

Xi. (“~ 




+ cos^Y^) ) 
^ ) 


V_ tan-l 

X ^ li c%. 


Therefore 


.splacemc 


4.65 


z(s) - z <*>1 + z: <t>i 

^•1 1*1 


2 °(i. 


* 2! ^ ® — ® 4* 


..Si ! '£ I ’T 


V i » 1 «4V 


C 2 ^ ' 

+ X { * “IP ^2 ^ ^3 I ^ X 

n ** 1 C ) X 


2 <S) 


ite AeceleratiOPt Z (S) 

2 « 2 «lS ^ 4>i= 

2 '15,— ip.e” + Z M,V 

1.1 1-1 * 
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(Y S ^ 

{ Of.v 6 # Jty, . Tia 2 \ 

0)0 4.^ ® V w { «4?a3oV i;; + U)o% ) 


where 

I- « 
4 


J 2. 2 XJs 

)2+ X 


X. sin<^ -^i I 


),' e “i® Sin -r, } 

+ < V )2 V i * 

I- I^IjOos (igaiLS ^ ^ ^ 

X '"V (I* 


5 , 


V 

C 2n TT ) ^ ) 

( j9-*)2+ (-ir"}2 ) 


(<^j2 «iS X 

, e -r cos X 

< V ^ 1 - 


iii) Contact force (F(s) 

P(s) » M ZiS) +ra Y(s> +(}44ia>g 

X 2 2 ^ 2 

-MX 2 %=l, Be->+2 


X i> » 1 


JL vl 


L_ (_2x <^\' 

i * 


— ^ I 4 4 cooig j . X 

3l*»l ( % ¥ 


( 2W1» | * 2 n^i^oos + <M+in)g 

‘ ‘ ) n.l 

4»2.2« gsaa^ai 

Following obsenrations may be maele CFig. 4*30 to 

4,32). 

1, In Fig, 4# 30# Fourier oonstants a^ and a^ have 
been plotted against n. Fixed pararaeters arelaa 
1* m 13si and h * Sam, It is seen that as ;r n increase s value 
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of &n decreased very sharply at upto n a 20, and then 
decreases slowly. Value of an at n a i is ^out 175 

tk 

times value of at n a 70. 

u 

2, Shape of the pulse is ctotained by two methods. 
First by directly taking assumed shape and secondly by 
considering equivalent Courier Series. A csWRparison 
between forms of two pulse shape obtained by two methods 
shox;s that two method gives identitltai iresults. 

3, On Fig, 4*32, ratio of absolute displ accent 
to the h has been obtained Jay two methods as in case of 
pulse shape, htere also, it is seen that the two method 
gives identical results. 

Thus it is concluded that response values chtalned 
by Fourier analysis are identical to that given by single 
pulse excitation analysis. 
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CHAPTER S 

RESPONSE OP PERIODICALLY OCCURRING RANDOM LOW JOINTS 
INTRODUCTION 

In the Chapter 4# we have ealcuXated all response 
paraineters deterministically by taking known of values of 
1 and h. In this chapter*, second order statisti<^ of 
absolute displacement* acceleration and contact fisrce are 
calculated by considering length 1 and depth h of the pulse 
as randcan* Such a analysis is useful in calculating the 
fatigue life of the structural cO!i^n<Kats« The occrurrances 

rji 

Of the pulses are not r andean* &ti4te -they occur periodically# 
Following assumptions are made# 

1) I<ength and depth are identically distributed and 
completely correlated random variables. They are related by 
the relation h » cl where c is a constant# This linear relation 
seems to be practicable as it is seen that increase in 1, 
increases depth h also (Pig# S«2) « 

ii> 1 Is uniformly distributed between 1^ anS 1^ as 
shown in fig. 5#1. 

Probability density function of 1 Is given by 

- V 1, 

Ijj • 1| * 2 


0 otherwise# 




i 

f i- 



FiGi S i 




I 


f 

h 






ii iz 

■ 


Ft 4 5'Z 
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® * 2 >1+ S ( usln(a(s-nL) + Y )*Q 

{ '*■ * A ( 

cos( SI (S—nL) +^) } ■*■ A ® ■**r ^(ucosY-^siny )) 

) I }) 

WHERE el « h 5 


Expected val«e of a function f<x) Is defined as 

eC£Cx»» ]^<H)fCx)ax 

where x is a random variable and 

(x) is the probability density function 
E(f^(x)) « jp3g(x)f^<x) d,X 

and variance of fCx) 

*^f(x) « E(f^Cx)) - |e f(x) f 


5.3 


5.4 


5.5 


using above definitions# expect®! value and variance are 
calculated. 

i) Expected value of ^ iS) t ECZ^^^CS)) 


Where Ptli is the probability density function hence 




i Vi 

1 . 


<S) 




dl 


17 ^ 


1 ^ 


5.6 


A closed fi^rm solution of the above integral is not 

p...ible a. i. a si.pXe functton 1. Thl. In^al 

is evaluated isimarlcally* 
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(ii) Variance of 

(S)> “ Y 1<®> 7 -*-, - 

J I2 - Ij 


1 ^ *“•' jL ^ 

2 1 




Hence, 


=^Pn,l(S) I - C )' 

1 < X ’ 

2 ( 2 J 

— ^ — f z* ,(s) ax - ( 1 r z„ .(s) at) 

I2- ( ij- li J ) 


i 

m. 1 ^ 

ifk •** «•' 


'n,l(s) 


"n,i(s) 


4 A 2 . 

z ® (s-iii.)+ |w :e 


1 S 1 


Cl 4^ . ^ 3 ^(S-nL+ « 5 t« 2 t 

— v{-aa)o« **■ X (Sl(s~nL)+Y) 

( X ¥ 

V 1.x 

-acos ( a (s-nL) + Y ) X + I u.* e ^.<s-tiL) + -)x ( a cosy 
X ) ^ 

-Q sinY ) X ) 5 

X) 

A) 2xpect,ed Value of f Z n,l^^^^ 

E{ Zq^j^CS)) ■ rVi: p “ = 

2 1 «j 


5.10 
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li) Variance of Z -(S)f (c^z* <s) 

UmI * 


E (Z*^- (S) ) « - ^ *2^ ^ 

^ n*l ' ^ ^ * I 2 - Ij ] \,1 ^ 


5.11 


hence 


n,i(s) 

JL 


"n,iis) 


# * i 


" 17 ^ 


« - ( 


< 1 


(ig - H 


• • 

^n,l (s) ) 


) 

di > 
5.12 


3. Contact Force» F_ . (S) 

Ilf 4> 

» M*Z^^^(S) ♦my* (S) + (mm)g 


5.13 


1> Expected value of F„ , (S) t E{F^ . (S) 

BfX Hfl 


E(F„ , (S)) « E(M Z„ , (S)) t E (» f • (S))+ (M4M)g 

JPif X Ilf X n 


*ME(Z^ ,(S)) t ME ( y„{S))+ (M4m)g 

Hf X H 

where E(S;^^(s)) = j^-j- ‘2„,i<s) ^ 

hence . 


TO 


X 


Xj- ^1 


J <SX t (M^)g 


5*14 
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ii) Varlaacs® of , (S) i ^ . 

«#1 n,l{s) 

J* 


r f \ ^{M4ia)g ) % 

2 1 . 

1. 

% 


lienee 

2 


2 

fct i 


^ ^n,l(s> • ® ^^UtlCsl^X 

Jtf% 

* *f‘ *l 3 [ f **** ^ CK4*Ei)g )%1 

2 2 j 


5.15 


ii(p||if villi’ ^||||e^ ^ 

rrr ( ‘"Vks ) ‘^® +{M-Hi^f) 61) 


5*16 


€b) ^ 


mt ilC» ;6- S ll£i •!• 


'b, 2iS) 


'lit2(S) 




^ V i 
ui ^ % 


< 1 **» © ^ I ♦ ^ B Cleogy -f u«ii»'f I ^ ^ %aa>»rm Qsinv % 

I 2 ^ A I K I 

(1 «© l[A»%tki0.is»m»}*'rl-¥ ^mn 

^a . AX 

^ Y I 1 + ^ • ^^®'****'^ |u8i»Y + aoo«r I X 5*17 

A i 

1) Sxi?«qfce6 Value of 


V3(S)* ■ i %S(S)® 


S.18 
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ii) Variance of Z (S) i ^ 

^ Z* t S J 

n ,2 

^2 

J <2<=> 


5.18 


hence 
X 

'n*2(S) 


a E \z\ AS)L Je (Z„ ^{S)) 


( p '“n,2 




f z2 

J n# 


< 

2 (S)- J ^ ^‘- 7 - 

di ( 2 1 


i 


X 


2 „^ 2 (S) di 


)2 

) 

) 

) 


5.20 


• # 

2) Absolute Acceleration! Z _(S) 
• n« 2 

2 

z ^(s) « Ct>^ o(^ 

l.*l 


n, 2 ' 


2 2 f 

^ cl : 2 u(s-nL) 

Z i V ( -utjOo e 


i SB 1 ^1 


) 


|~u.sin C^(&~nL} +Y ) «Q.cs>s( SlCS-nL) +y))+^U* 


) 


^U(S-lS|,) 4fl.C0SY ) ) 


5.21 


1) Eaqjecte f value of ^n, 2^^^^ 


1717 i V2(s) a,. 


5.22 


«* • 
f (S’ Z, 


ii) Variance of ° ^n, 2 <S) 

^2 

V2(s)> “ l^i l^ ' " J *'n.2(s) A 


5.23 
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henc» 


o •• 

Z. 


n, 2 (S) * ^ ^ 2 ^ ^ 2 (s)^ * ^ 


1^ \2 (S) - ci^rr^ 


t 

( 1 


2 ) 2 
•J ^ ll , 2 < S)^^j 


hen(» 


^ K,2( sP -“ irri :. 


*2 

i 


<3i $ 'M V' '^' 

2 - *1 


‘m,2(sr'^ i« - 2 


5.24 


3« %^<Sf,« ^n,2(S) 

'n.2(S) ■ "'^n.2(s) +»'’»»«S) *(H-m)g 

i) Esqpectsd value of ® 2(S) 


) 


®^^n,2(s)^ » E< M 2n,2(S) *** ® ^n(s) ^ ) 

« MeC * 2 ^^ 2 ( s )^ - HnE (* yjj ^ gj ) +( M^)g 


'2 .. 


where E( | 2n,2<S)^^ 


E( YJS)) » a *"“ ' f"' f 

« ^ 2 * n J 


2 *1 

i H 

^ Y ^( S ) dd 

n 


'1 

12 .. 


E( Yj}S)) «I ^ at 


5.25 


i *YL<S) at + (M»Wi5§ 


5.26 
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ii) Varl anc® of F^^^(s) t ^^n^2<S) 

® (^n.2(s)j- r: ii. S ''\,2(s)«i 




5.27 


hence. 


^n,2(S) * ^ <^n, 2 (S) ) -^0,2(8)^ X 

C ^2/ , ' I ^ ^ 

i ‘'«\.2(S) +”>"^n(S) + «««)9) )*' 

h X 1 j2 .. ' .. 

*" |7rpXj^ 1 <MZjj^2(s)4-inYn(S> + (M4«i)g><aiX 
Third Reaiont S xiL 4- 4 


5.28 


1) Absolut 


'n,3<S) 


2n.3(S) - Z„9(S-S„) + zsn h(S-S^) 


5.29 


Where g (S-S^) , Zg^ and h (S-S^^) are as defined in 
equations 4.25 4# 26. 4.27 and 4.28. 


y§1»g 9i ^0,3 (S) * ® ^^n.3 (S) ^ 

ig 

® C^n.3(S) )*“ i ^n,3(s)^'^ 


5.30 


\.3(8)* 


n.3(S)’ 


hence s \^ 3is) 

' M ^ #1 


rpq J 


2 


h,3(s) 


( ^n, 3 (S) I • C ® 3 (S^ J ^ 

^2 

< iji J %3(S) 

^1 


5,31 


5,32 



It 


2) Absolute Accel era tlona Z 


Ug 3 ^ S) 


^n,3(S)* -*• ^Sxfi 


5,33 


Mhere 


- ^0)0 


< 


g <s-s„) me V {s-s ) ( 


J (-2?U>ii^)+('?^U3|-COd^ ) 


(S-.S ) + <zJL 


sin ^ + (TiT^)^ <-2^:u3(f*3d)+Ut3tJo-a>d2j V 


) 


cos ~ (s-s ).} 

V o ) ^ aiyo^S^g ) ( 4 Clid, 2 2 2 

o SL_J ai JslnJ^{S-S ) (*$^031 0)4^) 


V* (S*-S ) a* 
h ' n' 


Tur 


- 2?oaJ% cos^ ) 


..lxpe.ct.ea .v alue....of 
® 3 (S) V * 1^^ ] 3 (S) ^ 


5.34 


# # 


li) Y;^jfB„an.3„ ,. g |, 2||,3(s)t 6*Z^ 


n,3(S) 


E ( Z.^ 


1 r2 *22 


‘n,3(s>^ “ j 


z-^3 (S) ai 


5.35 


hence 


‘^^n.SCS) * ^n ,3 

m 

‘2 4# 2 


*zf ,(s) > - 1® C2n,3(s>>J 


f - <1^ j’ Vsts)® ) 

JLm^ Am ^ ' ■* . ^ 


) 

^2 ♦# 


5.36 


2" n 






THESE ARE THE SECCiip. order statistics of response 
paraiaetera due to n^j^polse only^ Second order statistics 
for series of H pulses are obtained by Superinqposing second 
order statics of responds parameters duo to eadi pulse* Here 
also there id.ll be three regions* 
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(A) first Region! (N+l)L - | ^ S ^ (ItDL 


Absolu te Displacempintu Z(S) 

2(S) m y 2 4. 7 

nt% «'3<s) ^+l,l(s) 

(i) E3<pected value of 2<S) # £(Z(S)) 

N 

E(2(S) « 2: E(Z 4. K O 

nml «»3(S) T!i4l,l(S)' 

(li) Variance of Z(S)i g\{s) 

* 

N 

hence S-g^gj «E(2^(S))* lfi(2(S)) 


5.41 


S.42 


5,43 


5.44 


(2) Absolute Acdeleratlont Z^S) 


‘zCs) 


N 


"£ \,3iS} * ^+1,1 (S) 

B «» 1 

(1) Expected ^alue of *Z(S) i E C *Z<s)') 


N 


♦ ♦ 


E( Z(S)) a 2 E(z <S» 4. E( Zjj - j) 
Cli) Variance of Z(s)i ^2 ^(S) 


E(i*^(S)) 


H 2 

,^<'*n,3{s)> * 


> 


5.4S 


5.46 


5.47 
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ii€!iitoe s^2(s) ** ^ <K*(s (z{s)} ) ^ 

3) Contact Forcet F(S} 

FCs) m M 2*(s) 4-(M4ia)g "i* ra 

i) Expected value of F{S)* E(F(S)) 

ii:(F(S)) * ME(Z(S)41«E(Y* (S)) + (M4fli)g 

11) Variance Qf F(s)> 6%(s> 

E (P^ (S)W h\ (Z* * A ( '*Yj(S))+(M4ra) g) ^ 

+2tfe®(*2Cs) Y^<S) )+2MCM4«»)g * E 
(Z (S) ) t2m (M4nj) E ( * Yjj {S) )S 

hence dp(S) » E(F*(S» - (E<F<S3)? 

1 

(B) SECOND REGION! <N+l)Li: S£: (N+13D 4 ^ 

<l> Absolute Dlaplaceraent* z(S) 

N 

2i{S) - ^3(3)*** ®^S4U2<S) 

(i) Expected Value of ZC3)t E(Z(S) 

E(Z(S) « ^2 ^ "^"^^*1141,2(3)^ 

1 

■ 2 I 

Cii) Variance I of Z<S) 

" "H: : ' ■ ' ■ ' 2 ■' ' ' 

E{Z^CS)) » ^ ®^3»4l#2CS)^ 


5.48 


5.49 


5.50 


5.51 

5.52 


5.53 


5.54 


5#55 
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hence ^ 2 (s) “ ^ C2 ^(s))-.(e( z (s) ))^ 


C2) 

• ♦ 

2(S) 


Absolute Acceleration! z(s) 

^ • • • 

" ^■>.3tS)^-:Vl.2<S) 


DExpectea value of Z(S) i E( 2 *(S)^ 

N 

E(*Z<S)) « 2 E(*L „(S)) ^ I (1 ^(s)) 

n«l N+1,2 


ii) Variance of 2 (s) t (s) 

E(*Z (S)) « S B(*Zn,3(s)^ 
n»l 


^ ® ‘^n,3(s)*Vi.2(s)’ **“• 

m 2 ' 

hence G-^r^S) - E (Z (s)) - Ke (i(S» X ^ 


5*56 


5.57 


5.58 


5. 50 


5(»60 


(3) Contact Force t W is} 

P(S) » M Z(s) + Bt Yj^(S) 4* 5*61 

i) Expected value of F(S)# e(P(S)) 

E{PCS) « ME {Z(S) +WSC’Yj^<S) +(M4«)g 5.62 

li) Variance of P<S) I <5-p^(S) 

E(F^(S)) e M^(Z*^(S)) 4- «?£( Tf*(S)) 4* (M4ffl)g)^ 

4'2«ibB(Z* (S) t2MCM4aag EC’z(S)) 4-2»(«4*f!)iryjj(S)} 5,63 

hence ^F^{s) • bCf*(s)) * ({e (fCs>))^ 


5.64 
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(C) THIRD HSGIONf S»{N+1)L 4- 4 

2 

1) iUssolute Displacement! Z(S) 


Z(S) 


H+1 


n*sl 


z 


n,3 (s) 


5*65 


i) Expected value of Z(S) * £(Z(S} 


E(Z{S)) « 2 E(Z 

litti 


n»3(S) 


) 


±1) Variance of Z(s)t ff^Z(s) 

E(2:*(S»= X n.3(sp 

na*i 

hence e^z|s) » E{Z^(S)) - (E(Z (S)j^ 


5*66 


5.67 


5.68 


* • 


(2) Absolute Acceleration! Z{S) 


(S) 


H4*1# . 


:z %3(s) 

i) Expected value of i E(Z 

N+i 

E{z*{s))« Y ® {z*n,3(s)) 

ii) Variance of Z* (s) i 6^^ (s) 

N+l 


« # 


E('Z^{S)) » Z S ( Z 

IS«1 


I^(S) 


) 


hence 6Z^(S) «. E(’zfs» - {E< Z 


5.69 


5,70 

5.71 


3) £».te.fe..la£gfii. 

F(S) « M*Z{S) +CM4®)g 


5.72 
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1) Expected value of P(S)i E(F(s)) 

E(FCS» « METziS)) + (M+m)g 
ii) Variance of F(SJ*» 6^F{s) 

E(F^(S)) « M^E(*z2(s); +2M(M-Hn)g E(*2(S)) 

hence (S^F(S)«. ECF^(S))~ <E(FCS)))2 

s,2 isayi.tra 

Results obtained in this section are for vehicle 
described in Tsfol© 3*1, other data used are- 

i) Maximum length of pulse » 7.5m. 

11) Minimum length of pulse (1^) « 1.5m. 
ill) Mlnlinum depth of pulse « 5ram 
iv) Maximum depth of puise ** 25OTa 

C » 1/300. 

Following observations be mede (figs. 5,3 to 5.8) 

1, Expected value of absolute displacero^t has been 
plotted against S/1 for different velocities in Fig. 5. 3. It 
is seen that expected values are high for low velocligf. At; a 
velocity of 30 kHv'hr. Maximum value is 25 mm while it is only 
S.Smm for a velocity of 15€ac3%%r» 

2« R.m.s, value of absolute displacement has been 
plotted against S/1 for different velocities In Fig, 5.4. 
There also r.m.s. value are high for low veiicity. Maximum 


5.73 

5.74 

5.75 
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value of r.m.s, value at velocity of 30 km/hr is I 4 n 3 m 
while it is 3«0 mm for a velocity of I50kro/hr, 

3. Ratios of expected value and r.m.s. value of 

absolute acceleration to the acceleration due to gravity 

have been plotted against S/1 for different velocities in 

Figs. 5,5 and 5.6. It is seen that the ratios are high for 

high velocities. At a velocity of 30 km/hr, max value of 
# ♦ 

E( Z) is about 0.06 times g and 62 is 0,09 times g. For a 
velocity of 150 km/hr, corresponding values are 0.14 and 
0.17. 

4. In Figs. 5.7 and 5.8, ratios of expected value 
and r.m.s, value of contact force to the total weight of the 
vehicle have been plotted against S/i for different 
velocities. It is seen that the increased velocity higher 
maximum value of the ratios. At a velocity of 30 kn/hr, 
maximum values of E(P) and ^F are 1,07 and 12 times the 
total weight of vehicle. For a velocity of 150 kny'hr, 
corresponding values are 1.09 and 0,20 respectively. 

Thus it is concltttlcc*. that the incrdase velocity 
lowers maximum r.m.s, and expected value of absolute 
displacement, while opposite nature is seen in the case of 
acceleration and contact force i.e, increased velocity 

increases maximum r.m.s. and expected values, 

♦ 
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CHAPTER 6 

RESPC^jjiE OP VEHICIiE TO COfSTINUOUS AND TOTAL RAIL TRAQC 

UNEVENNESS 

INTRO-DUCTICXm 

The continuous unevenness of the rail trade can 
be treated as hcwio^nous random process in space coardlnates 
when a vehicle moves on such a track with constant velocity, 
the base excitation is stationary. The process is asstamed 
to have zero mean and spectral density given by 

$yy(0-) - iQUae 

IT + 

» 0 . 

hhere Q. is spatial frequency, c and f are constants* The 
roughness constant C is given by 


where S is the r.ra.s. value of ground unevenness, 

Secemd order statistics of response parameters are 
calculated in section 6«l4 In the section 6.2, continuous 
unevenness is superimposed upon descrete unevenness due to 
random low joints. Second order statistics for total 
uneveiwess is obtained by assuming that two process are 


uncorrelated. 
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6.1 Oa^DER statistics Of 

TmCK UNEVENNESS 


RESPa<.SE TO CONTINUOUS RAIL 


Equation of Motion* Equation of motion for single 
degree of freedom model is 


+2^U)oZ +Ud| 2 « 2TgU)oY + U3oY 6.1 

for oonstant velocity motion S » Vt» above equation becomes, 

V^r.« + 2'^Vu3oZ* » 2?U3 oVY' + u)2y 6.2 

Solutioii Of equation of nration* 

The solution of equation 6,2 be expressed as 


2 

2(S) * 

oo 

f .TO o 

A e V + 1 dF(Q.)e^^^^H(a} 

6.3 

i«»l 

-00 


where 

o<^ » ( 

i « 1,2 


H(S1) 

•» (1^ IkJoJVQ^ + U)h 

614 


(a)|-v^^+2J^avu3o) 


P are constants and to be deteriained frcan initial 

JL 

conditions. 

■oO 

y<s) » J e dPCQ) 

^00 

and PC^) 1» related to the P.S.d function ofY(S) by 
the folloving equations 


E (IdP(Q)i^) *fyy{QH^ 

6.6 

also E(dF(Q» » 0 

6.7 

as E (Y(S3} « 0 

6,8 

for aero initial conditions, we can write as 


oo ^ 

Z(S) « J 

6,9 

and*&cs) • - J df(a)H(Q) 

6.10 

^OO 

Jf 
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1 ) iixpected value of Z(3) aad*2{S>i E(Z{S)), E<*z(S» 
00 

f JQ.S 

Eiz)(s))m J H{QL)e E(dP(a»»0 

-00 

frcMit eq ation 6,7# and 

E( Z{S)) w (Z{S)) m 0 from above. 

dS 


C2) Qavariance of Z(S) and Z(s)t 


•ihG response and its defivaties can be closed as 
an cscillcxy process. Hence the evolutionary spectral density# 
which is the common spectral density in present case can be 
determined by the use of following esqsression given by 
Priestley. 


^zzJLQ.) m ^yy(Sl) HCQ) H*(a> 6,11 

Where i*) denotes conjugate operator, 

4 s the eigen values and eigen vectors occ^r in coR^lex 
conjugace# alcove expression can be written in a sintpler 

form. 


^ZZ{Si.) m ^yj/CQ. )H(Q.)H(-S1) 6,12 

covariance funi^ion for the response is 

K (sr, S2) . 1 

with the same argument, the conjugate operator may be 

dropped with chatfig® i” * Kence 

tS fO 6.13 

K^^CSl, S2)* J$yyCCI)H(a>H(-^ J e 
SB*' «« 

* uiO.\ HC-C2) and^yyiO-) we get. 
Substituting value* of ■ 
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K(sl,S 2 } = 
Z!S 


J 



(2? U)o yjSI » coj) C-2g(A3ovjQ + co j} } 
C U)o« ^^‘f2?coo ( <lOo*V^ 


• ® 6*14 

111© Impx'c^eir Integxal aliove can. b© ©vain ©ted fey 
contour int^ogration with residue theorem* There are two 
polos in the upper pjirt of coittplex domain at a. Jr and 

and m i « 2 

cs,^ H +\/ ^ ^• 1 ) U)o 


both the poles are of first ord,er* Hence 


K„ 

za * 


f 

C 2 TT J{ 




2 2 2 2' 
41 ? ‘V a)oT e 

"■ ■■■'’■ "' ■ ■ ■ m ' "A »»«■ '■ " "I" 
.2 .»» 2„2 


-y(si-s 


2 ) ) 


( C 2dr{ cbo+v-T - 2?aJoW) ( (i)Uv 




( 


t 2 i 5 u 3 ovr)) 

( colt (aj/'S^-i) 0)0 V 

^ 0(^1 ( u) 0 * 4 ' 5 ^cx)o^ °^ 2 ^) I 

* !!|! ) c a)(?+ «2 - 2^ a)o«2) ( - 2J/^-l) V 1 


« "' 1 C 


I cuA- 

< (X(6j 3§+V^^'^)^- 4 -S^COo^V^r^)) 
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Above equation Inaicates that lc^^(Si,s2) is a function 
of Sl-;i2. Hence displacement and its, time derivatives 
are homogenous process* In that case,. 


also If Sl~S2 ** S then 


R, Cs) a nc X A y?)....6. 

^ ir ico?+v • 4 


« C COq- 4 ^^0)1 


Cic^~ •■•i'y-) i (00+ 2'S^O®^!} ( 

^ e {(0^- 4 “S^tOo^oc®) X 

«! 2 ^ 

) ( COo+% * 21?U)oof«) C 

t m- W 

¥ 

■Iharefor. ^ Ico^^^VcOoV; 

.<.2<s,) . K^^(o) - ac I 

* (U?- 4 -i’uo’ «, > 

,< i i ii ,i . i.ww« i *i > <> > i y .ww jiifl.:Cii[f»ii.iirf 

2i-) ca)o+°'i - 
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. ( 4 Li)n^ ■ \ 

) ( mI+0(|. 2-5i0o«2) ( J 

2 

Now variance (S (s) is 

^2^(&) * E( 2^(8)) * ((E (2 (S)))^ 

3' 

m 6» 17j as E(zCs) is zero* 

Since Kjgjg(Sl,S2) is a function of (Si-SZ), then for 
evaluation of the covariance function of the second 
derivatives# 'che differentiation has to be t^:en with 
respect to (S1-S2) « Hence 


K'* *» C£*l#b2) 

xz 


(S1-.S2> 

^•£S*(S) where S a ei-.S2 


^ I ( 2 ^ U3oVil J.» col ? 

-flo ^ ( U)|- V%l®t2?u)oVJ-SI) <<jOo-V^^-2tS vacdop.) ) 


dP 


wChsreforc# E{ 2 Cs)'} * Co5 

Hence variance ^2 Cs)i 
(S**2 


6.18 

6.19 


Hb) 


I( I ^(3))- X3(Z (3))X 


.1 e|*z2(s) Xas s CZ<S))is sero. 


6,20 


• • 


3) Ctovariaaace of Y(S) aiKi 3 
w# have ^ 

CO 

K CS1.S2) « J $yyCQ.) 


-00 


Cs )* 5 ( s}s 
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jSl{si«S2)da 


-CO or 4- T 


Vy* 


But K. • 

yy 


Thcrcfc 


and EC 


Thc^refc 

6"* • 2 ^ J 


E C 


Again %» 


Z(S ) . 

A 


Y(S^) • 

m 


6* 21. 


also K..^*{ 8 U. 2 ) * as , 2 ^ s ,2 


but sin 


henea 
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Since is also a function of (S3-.S2) therefore 

Kr'y(Sl,S2) = Rj-f (S1-S2) » K2^(S1.S2) 

^ ( 25 u)rf^SlJH.a)l)e aa 


i 


qq: 


• Q-o 


dT? (aio-¥^il^4-2?o)oJvSl ) 


hence- R^* y* ^gj 


( 2XyoVJ^l dSl 


6.27 


(*) gxpeeigt j ,,al a ti_ag A -ya£A.§a£e....^^ ,fec^^Cs)t 

F(S) as H Z(S} 4- ra Y (S) + (M4iR)g 6.28 

hence 

1) Eatpec'ced value of P(S) *» E(F(S))i 

E(I (S)) a ME(*ZCS)) + WE (Y*(S)) + (M4ro)g « {M-iTO)g 6.29 

as E (Z*(S))and E(*Y(S)) are zero. 

ii) Varlanca of P(S> t ^P (S) 

E{P^(S» « M^E(*Z^{S)) * ro^E ( Y^<s)) + (M+ib) ^g^ 

4.2 Mi!iEC* 2CS) *YCS)) +2CM-m)J^ E( S^S)) +2(M4iB)nigS { Y (S)) 

Si M^E(*2^(S)J + inmd g^+a«n!B( Z(s) yCs)) 6.30 

by substituting valuea of £ (z’^(S)>,E Xfhs)) and 
E( Z* {S)*Y(S)) from above. E (S' ^(S)) can be calculated, 
hence. ®’P*(S) * E( P^(s)) - XJ3(P(S))X 
•W^C*2^(S)) tll^l(*Y^<S)) +2M»B( Z(S)*f (S)) 


6,30 
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6*2 SECOND OivDER STATISTICS OF RSSPOSSE PARAMSTERS TO 
TOO[*iVL Ui3;iVENNESS| 

In the last section# we have obtained second order 
statls-tlcs due to continuous unevenness only. In this 
section# continuous unevenness is superinaposed upon descrete 
unevenness due to low Joints, Second order statistics of 
response process due to total unevenness is, obtained by 
assuming the two processes are uncorrelated. 

Let Yj^(s) be the process due to continuous unevenness 
of rail track \ tr uted as homogenous process with zero mean, 

Z^Cs) the response process due to Yj^(s) and can also 
be treated as htMBogenous process with zero mean, 

Y (S) be the process due to discrete unevenness of 
random rail joints, 

Z (S) be the response process due to Y^Cs) 

By superpositioning on principle 
Zi3) li Z^(S} + Z^ts) 
and Y(s) m Yj^(S) tYjCs) 

1) Second order Statistics of ZCS)i 
i) Ej^ected value of 2 |S) i E(z(s) 

E(2)CS))fc E(Z^(S)) + E(Z^{S)) « 
as E(2j^(S)) m 0* 

ii) Variance of 

BCZ^CS)) - (B(zJ(S)) t(E(z|(s))4E<Z3^(S)Z2CS)) 6,32 
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6.33 


E (Z^(S) Z^CS)) a 0 ae the two processes have been assumed to 
be uncorrelated 


^Z^(S) «ECZ^(S))- XE(Z(S))X^ 

2 

m e{z^(s)) +e(^«(s)) - 

JL W mi: 


as 


E(zJ(S)) + S^ZjCs) 


where ®Z?{S) » E(Z^(S)) — X^(^a(^))X 


6.34 

6.35 


2 ) Second Ord er StaUstlcS 0 .f_ 2 (S) 


*2iCs) ss 2 jC5)+ ZjjCs) 


i) Expected value of Z{S) i E{Z(s)) 


E( Z(S)) « E( Z1{S)) + E( 

# • « * 
a SC as E( 

il) Va lance of zis) t6Z ^{S) 

ji{*i*{S)) - -(‘zJ(S)) t E('2a(S)) +2E(Z^(S) Z^fS)) 

. E{’Zj^*(£))+ Zj(S)) 

heno,ffP(S) .e{*£»(S))-XE(‘z(s)) X* 


.E('z*(s))+E(*z*ts))- xzc’vs)) X^ 

a E( Z^cs)) 

where - EC'z^Csn- «S (’^zts)) X 


6,36 


6.37 


6,38 
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{%) Second Order Statistics of Contact ^orce P(S) 

P(5) « m*2{s) 4- m*y(S) +(M4in)g 

1) EsQiected value of F{S)* E(P(s)) 

E(P{S)) = I-1EtS(S))+-«rtE(*-Sr(S)) +{M-Hn)g 6,39 

where E tS(s))as calculated above, and 

E(*Y(S)) a EC^^iS)) +E Vy^Cs)) » E(Y*(S)) as 

E{*Y3^(S» » 0 6.40 


(ii) Variance of ]f(S)j 6^F(S) 

E{F^(S) « M^E(*Z^(S)) 4ffl^E{*Y^(S)) +(M+ni)^g^ 


•f2M(M-«n)g E(2(s)) +2irt® ( 2(S) Y<S)) + 

• ♦ 

E( Y(s)) 

where 


(2*^(8)) « E(*Z?(S)) +Ei 2,(3)) 


E(*P(£)) « 2{*y|(0)) 41 2E( Y^(S) Y^CS)) 

« E(*Y?{S)) + 


6,41 


6.42 


6 ^ A3 


ajB*Y-CS) and^Y-iS) are unoorrelated process and therefore 
* * ' • # 

E( YlCS) Y^iS)) m 0 

E C ZC^)) m E( 22(6)^ 

S (YiS)) » ^ ( Y2{s)J 


6,45 
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E(*Z(S)'y(S))« E{Zi(S)\cS)+l^(S)‘Y2(S) ;J2(S) 


Y- (3) *}- Z (S)*Yo(S)> 


« E(*23L<S)*i-i(S)) + E(*Z2(S) \i S)) 
hence (3) a E( F^(S))- %EiF{s)yj^ 

6*3*Hi-:SUhTS 


6.46 

6.47 


Resul'ts obtained in this section axe fox a vehicle# 
the data ,^.or i.hich is given in Table 3,1* In aaditicn to 
that other data used are 


Danp lng coefficient 0.2 
Roughness constant C e 0.000Q5 
T w 0 • 1* 

Following observations may be made (Figs. 6.3 to 6,8> 
1, In Fig, 6.3# ratio of the r.in.s, value of ^solute 
displacement to the r.in.s. value of ground imevenness has 
been plotted against velocity. It is seen that the ratio 
increases linearly with velocity being alaout 1.2 times the 
r.m.s. value of Y at a velocity of 150hi^/hr, 


2. In Fig. 6*4# ratio of the r.ra.s. value of absolute 

acceleration to the acceleraticn due to gravity g# has been 

plotted against velocity. At lotv' velocity# ratio increases 

sharply# but for higher velocity ( > 601cm/hr) • Variation is 

# • 

approocimatoly linear. At a velocity of ISOkm/httCZ is 
0, 27 times acceleration due to gravity g* ■ 
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3) In Flg« 6*5 ratio of the r*m*s* value of contact 
force to the total weight of vehicle has been plotted against 
velocity ratio increases very sharply at low velocity* At a 
velocity of l50kro/hr,SF is 0.24 times the total weight of 
vehicle. 

4) Assuming the unevenness to be normally distributed 
the 36 values which have the excedence prcb^lllty of 0*0027 
are 36^# 36’*^ and 36^. 3d values have been plotted in Figs. 

6i6* 6.7 and 6.8 against velocity. 

5) *vhen discrete unevenness due to low joints is 
superimposed upon continuous unevenness (section 6,2)* the 
second order statistics of the resulting process is obtained 
by assuming the tvK> processes are uncorrelated. Since the 
expected valtaes of displacement and acceleraticxi process due 
to continuous unevenness are zero* the expected values of 
displacem' nt and acceleration process due to total unevenness are 
same as that due to di^screte unevenness. R*ra*s. value of the 
total response process is obtained by using following relati<xi. 

0" « (dj + 

where <y. and<y, are the r.m.s. value of response process due 
to descrete and contlnuoxis unevenness. Expected value of 
contact fot'ce due to continuous unevenness is the weight of 
vehicle. R.m.s. value of total response process has been 
plotted in Figs. 6.9* 6.10 and 6.11* against S/i for velocities 
of 30 and 60 mnAir. It Is seen that increase velo«dty gives 


iligher values^ 



Thus it is concluded that in the continuous 
unevenness# r*m*s« values of acceleration and contact 
force increases with velocity* Results very iiuch depend 
upon cho en value of roughness constant C. Expected value 
of total response process is saoae as that oi response 
process due to discrete unevenness* R*itt*s* value of total 
response process increases with velocity*' 
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CHAPTER 7 


SUMMARY AND CONCLUSIONS 


7.1 SUMMARY 

A vehicle moving in contact with the gro.tnd follows 
the ground profile. The ground profile# therefore acts as a 
base excitation for the vehicle through the ipoints of contact. 
Present investigation deals with the response of railway 
wagons to the vertical unevenness of the rail track. Two 
types of unevenness of the tra<^ have been considered. First 
continuous unevenness due to the irregularities of the rail 
profile and second, \anevenness due to hogging of the rails 
at the joints* The vibration effects of each type of the 
unevenness have been analysed separately and superimposed. 

Idealised models have been constructed for veMcle 
and track unevenness. It h^ been assumed that vehicle 
moves with constant velocity. The governing equations of 
motion have been decoupled to give set of first order 
differential equations. 

Response statistics for periodically occurring low 
joint have been obtained by (l) single pulse excitation 
method and (il) Courier analysis and con 5 >ared. 

Second order statistics of displacement# 
acceleration and contact force for periodically occurring 
random low joint have been obtained for fatigue analysis# 



Second order statistics of response paranieters 
have been calculated for continual s unevenness of the 
rail profile* ITie continuous unevenness has been treated 
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iv) (a) for a velocity less than 50 knv^hr* 
interaction effect of only previous pulse is significant, 

(b) jPor a velocity lies between 50 to 100 kni/hr, 
effect of at least 3 pulses should be considered, 

(c) S'or a velocity greater than lOO ki^/hr, 
effect of at least 5 pulses is necessary in smaraation. 

This interaction effect villl again be insignificant as 5 

rail lengths are vjelded for hl^ velocity as found in practice, 

v) For a design DLF of 1,77 and at a velocity of 60JcinAir# 
probability of failure of the joint for one 3cro length of tracSc 
is as high, as 0*96, It increases as we decrease the value 

of design DbF, 

(B) Response of vehicle to continucus unevenness shows 

that- 

i) r*m*s, value of acceleration and contact force 

increases linearly at hi^ velocity. Results very 
BRidhi depend upon chosen value of roughness constant, 

(C) Response of vehicle to total \mevenness shows that- 
1) Expected value of total response process Is same 

as that of esqpected value of response process due 
to 4*screte unevenness, 

ii) r*m,s* value of the total response process increases 
with velocity# 
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7,3 FURTHER WORK 

The work can be extended as followsi 

^ A) Vehicle llodel i 

i) Single degree of freedom model can be replaced 
by degrees of freedom model vdth hawe and 
pitch degrees of freedom, 

ii) Interaction effect of successive passage of the 
two axles of vehicle can be considered, 

ill) Parametric study of natural frequency (W)o ) and 
dan^.l-ng ratio (?) can also be considered. Non 
linearity in the spring and dsrsper can also be 
taken into account, 

m irack Model » 

iv) low Joints have been assumed to have frozen 
profile, The redls can be assumed to be continuously supported 
on flexible foundation with the weic^it o£ the ties uniformly 
distributed and added to that of rail, 

v) Length and depth of the joint (l#h) have been 
assumed to be uniformly distributed end linearly correlated, 
■By collecting statistical data, actual prob^ility 
distribution can be determined, 

(C) Vehicle Motion ; 

vi) Constant velocity raoticsi can be replaced by 
variable velocity run, bue to change in space time 
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relationship, the governing equation of motions in space 
coordinates will change* 


vii) A detailed fa'tigue analysis of tlis vehicle and 
trade coirponents can be considered* 

viii) Reliability based design of the traeik can be 


done 
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APPEHDIX^I 


REDUCTIQEJ OF EQU^Xm 2. 2 A SYSTEM OF FIRST ORDER UKCOUPEED 

SQUATIQKS 

Consider equation 2,2 let 

hi = \ l I = {a-silY +u)|Y ) ^-1 

iO 1 > {-1 0 ) 

M a ( . ) , K » ( , ) A-2 

(1 2’SU)o) <0 U)^ ) 


Equation 2,2 can now l>e written as 

(M) tqt* (K) XqX« (i(t)) A-3 


The solution of homogenous part of equation is assumed to be 


of the form 



Xq X*^ e A- 4- 

and (D)(4)) 

= 

Where 



D s - » 

C -1 0 ) -1 

<0 u),?) 

< 0 1 

<G 1 > 

^ 1 2?U)6 ) 

) 


) 

) 


A-5 


Ai-6 


Where « and 4> are the eigen values and eigen vectors of D, 


The model column matrix is 

'2 


(§) 


* I 

XC4> ) X 


Ct> 


2 

t’2 


A-7 
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l.et Iq Ct) X« t ^ t XXCtjX 

Substituting into A»3 and using orthogonal property 
of eigen vector* yields a system of first order equations* 

tat i ** ^ ^^«»9 

Where Ca^CdoY + u}oV ) 

m 2 C °<j. t ■§ ti3o > 



